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DESORPTION  AND  SURFACE  DIFFUSION: 


NITROGEN  ON  TUNGSTEN  (110) 

Anthony  J.  Polak 

Coordinated  Science  Laboratory  and 
Department  of  Metallurgy  and  Mining  Engineering 
University  of  Illinois  at  Urbana-Champaign,  1977 

ABSTRACT 

Scanning  Auger  spectroscopy  has  been  used  to  monitor  the  migration 
of  nitrogen  atoms  down  a concentration  gradient,  on  the  (110)  plane  of 
tungsten.  This  concentration  gradient  is  formed  by  the  electron  stimu- 
lated dissociation  of  molecular  nitrogen,  initially  adsorbed  at  95°K. 

These  deposits  are  found  to  be  thermally  stable  to  several  hundred  K and 
uneffected  by  further  electron  irradiation.  The  adsorbed  atoms  form  two 
different  surface  structures,  depending  upon  the  surface  concentration. 

The  desorption  kinetics  of  both  molecular  and  atomically  bound  nitro- 
gen have  been  studied  by  a combination  of  Auger  spectroscopy  and  flash 
desorption.  The  desorption  of  atomically  bound  nitrogen  occurs  as  a second 
order  process  with  an  activation  energy  more  than  7 times  the  activation 
energy  of  diffusion.  Despite  the  large  binding  energy  of  the  atomic 
nitrogen  to  the  (110)  surface,  153  kcal/mole,  there  is  virtually  no  change 
in  the  work  function  upon  adsorption.  In  addition  the  saturation  coverage 
of  the  irradiated  atomic  nitrogen  is  found  to  be  more  than  twice  that  of 
the  nitrogen  adsorbed  at  room  temperature.  Both  are  found  to  be  in 


identical  binding  states. 
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I.  INTRODUCTION  AND  HISTORICAL  REVIEW 

The  stability  of  chemisorbed  layers  is  of  interest  in  understanding  the 
chemical  processes  that  occur  on  surfaces.  Such  technologically  important 
reactions  as  catalysis,  sintering,  grain  boundary  segregation,  and  thin 
film  formation  depend  critically  on  desorption  and  diffusion,  the  two 
mechanisms  governing  the  stability  of  chemisorbed  species.  Not  only  is 
diffusion  important  in  surface  reactions,  but  it  is  also  of  considerable 
theoretical  interest;  the  magnitude  of  the  barrier  to  atomic  motion  provides 
important  insights  into  the  nature  of  the  gas-metal  bond.  With  this  in  mind, 
we  have  studied  the  diffusion  and  desorption  of  a gas  chemisorbed  onto  a 
single  crystal  plane,  specifically  nitrogen  on  the  (110)  plane  of  tungsten. 

Although  the  diffusion  of  metal  atoms  on  single  crystal  planes  has 
recently  been  opened  up  in  investigations  using  the  field  ion  microscope  [l] , 
relatively  little  is  known  about  the  diffusion  of  gases  chemisorbed  on 
metals.  Comer  and  coworkers  [2,3],  using  field  emission  microscopy,  studied 
the  diffusion  of  a number  of  gases  on  the  high  index  planes  of  a tungsten 
field  emitter.  They  deposited  the  gas  on  one  side  of  the  field  emitter  and 
upon  warming  observed  the  spreading  of  the  gas  layer  over  the  surface.  From 
the  rate  at  which  the  gas  covered  the  surface  an  estimate  of  the  diffusion 
coefficient  was  made.  For  oxygen  on  a tungsten  field  emitter  the  barrier 
to  surface  migration  was  found  to  be  30  + 1.5  kcal/mole;  for  hydrogen  on 
tungsten  this  barrier  was  estimated  at  9.6  - 16+3  kcal/mole.  The  activa- 
tion energy  of  diffusion  for  hydrogen  on  tungsten  amounts  to  approximately 
207,  of  the  binding  energy,  depending  on  the  nature  of  the  sites  involved. 
Using  the  same  technique,  Ehrlich  and  Hudda  [4]  followed  the  migration  of 
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^ nitrogen  over  the  high  index  planes  of  a tungsten  tip.  From  the  spreading 

I out  of  the  nitrogen  layer  they  estimated  that  the  barrier  to  diffusion 

was  ss  20  kcal/mole  around  the  [ill]  , and  ss  35  kcal/mole  for  the  sharply 
I stepped  region  surrounding  the  [lOO]. 

Very  little  work  has  been  done  on  single  crystal  planes,  and  only 
^ recently  have  quantitative  studies  become  available.  After  our  research 

j was  underway,  Butz  and  Wagner  [5]  reported  the  observation  of  oxygen 

diffusing  on  a macroscopic  surface.  They  used  a masking  technique  to 
create  a deposit  of  oxygen  on  the  (110)  plane  of  tungsten,  and  contact 
potential  measurements  to  follow  the  migration  of  the  gas  on  the  surface. 

With  this  technique  they  found  a concentration  dependent  diffusion  coefficient 
which  can  be  described  by  a constant  activation  energy  of  27  + 2 kcal/mole 
in  the  coverage  range  0.4  to  0.9  monolayer.  Wells  and  King  [6]  , on  the  other 
hand,  using  a molecular  beam  nozzle  source  to  deposit  oxygen  on  the  (100) 
plane  of  tungsten  and  secondary  electron  emission  to  detect  the  oxygen  on 
the  surface,  were  unable  to  observe  any  diffusion  of  oxygen  along  the  surface. 
These  are  the  only  studies  that  have  dealt  with  the  diffusion  of  gas  atoms  on 
macroscopic  surfaces. 

The  overriding  obstacles  in  the  study  of  gas  diffusion  on  macroscopic 
surfaces  are  1)  the  development  of  a technique  with  sufficient  spatial 
resolution  and  sensitivity  to  detect  the  variation  in  surface  concentration 
that  occurs  as  a gas  diffuses  along  the  surface,  and  2)  the  actual  creation 

of  a surface  concentration  gradient. 

f 

[ 


/ 


3 


These  two  problems  can  be  overcome  by  the  use  of  Auger  spectroscopy, 
as  was  done  by  Jannsen  [7].  Using  electron  stimulated  desorption  to  create 
the  localized  deposits  and  Auger  spectroscopy,  he  studied  the  diffusion  of 
sodium  on  MgO,  for  which  the  barrier  to  surface  diffusion  was  found  to  be 
< 1.2  kcal/mole.  We  chose  scanning  Auger  spectroscopy  as  the  technique  by 
which  to  study  the  variation  in  surface  concentration  of  an  adsorbed  gas 
layer;  Auger  spectroscopy  provides  us  with  an  extremely  sensitive  surface 
tool,  with  potentially  excellent  spatial  resolution.  The  Auger  electron 
gun  can  also  be  used  to  create  the  localized  surface  gas  deposits. 

It  is  well  known  that  the  characteristics  of  the  gas  adlayer  will  depend 
critically  upon  the  exact  gas  and  crystal  orientation  being  studied.  There- 
fore the  use  of  a single  crystal  cut  to  expose  only  one  plane  is  imperative. 
The  choice  of  the  crystal  substrate  is  restricted  to  elements  which  have 
minimal  surface  impurities  and  surface  defects,  can  be  cleaned  in  an 
ultra  high  vacuum  and  are  able  to  withstand  the  repeated  bakeouts  necessary 
to  achieve  an  ultra  high  vacuum  environment.  The  only  materials  that  meet 
these  requirements  are  the  transition  elements,  which  in  general  are  easy 
to  prepare,  clean  and  observe  in  single  crystal  form. 

Tungsten  was  chosen  as  the  substrate  material  because  of  the  extensive 
research  which  has  been  performed  on  it  with  a number  of  adsorbed  gases 
[8-12] . Nitrogen  was  selected  as  adsorbate  since  it  is  the  simplest  diatomic 
gas  that  can  be  detected  by  Auger  spectroscopy  and  because  of  its  unusual 
behavior  on  the  (110)  plane  of  tungsten  [13,14], 
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The  behavior  of  nitrogen  on  the  (110)  of  tungsten  is  quite  different 
from  that  on  the  ether  planes.  The  most  striking  differences  are  in  the  rate 
of  adsorption  and  the  work  function  change  that  occur  after  nitrogen  is 
adsorbed  on  the  surface.  Typically,  the  initial  sticking  coefficient  for 
nitrogen  on  planes  other  than  the  (110)  amounts  toss  .33  at  298°K,  rising 
to  sa  .51  at  115°K  [8] . On  the  (110)  plane,  however,  the  sticking  coefficient 
is  only  . 004  [ 13]  . 

The  variation  in  the  work  function  of  single  crystal  planes  with 
different  surface  coverages  has  been  studied  by  numerous  techniques,  on 
both  microscopic  and  macroscopic  samples.  On  the  (110)  plane  of  tungsten, 
Delchar  and  Ehrlich  [14],  employing  the  Kelvin  contact  potential  method,  were 

unable  to  detect  any  change  in  the  work  function  after  room  temperature 

-2 

exposure  to  nitrogen.  Only  by  increasing  the  pressure  to  10  Torr  or  lower- 
ing the  temperature  to  190°K  did  they  observe  a change,  which  amounted  to 
-.13  eV.  From  the  pressure  dependence  of  the  work  function  at  room  tempera- 
ture they  deduced  that  the  nitrogen  was  bound  with  < 9 kcal/mole.  Using 
field  ion  microscopy,  Ehrlich  and  Hudda  [4]  were  also  unable  to  directly 
observe  adsorption  of  nitrogen  on  the  (110)  plane.  Even  at  high  exposures, 
when  changes  were  noticed  on  all  the  surrounding  planes,  adsorption  of  nitrogen 
on  the  (110)  plane  was  not  noted.  Only  after  the  emitter  temperature  was 
lowered  below  190*^K  was  nitrogen  adsorbed  on  the  (110)  plane.  This  result, 
however,  should  be  viewed  with  caution.  The  inability  to  obser'/e  nitrogen 
may  have  resulted  from  the  desorption  of  nitrogen  adatoms  from  the  emitter 
by  electron  tunneling  from  the  imaging  gas. 
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In  LEED  studies  by  Stern  [15]  and  by  Dooley  and  Haas  [l6l  no  change  of 
the  clean  tungsten  (110)  pattern  was  found  after  exposure  to  nitrogen. 

These  results  were  interpreted  to  mean  that  nitrogen  did  not  adsorb  on  the 
I (110)  plane  of  tungsten  at  room  temperature. 

The  desorption  behavior  of  nitrogen  from  the  (110)  plane  of  tungsten  has 
been  found  to  be  similar  to  that  from  the  other  planes  of  tungsten.  At 
I temperatures  less  than  120°K,  nitrogen  is  adsorbed  as  a molecule  in  the  y 

state.  Nitrogen  bound  in  this  state  desorbs  with  first  order  desorption 

* kinetics  at  temperatures  greater  than  150°K.  At  higher  adsorption  tempera- 
tures, nitrogen  is  bound  in  an  atomic  state,  referred  to  as  the  5 state. 

From  the  (110)  plane,  the  3 state  desorbs  following  second  order  desorption 
kinetics  with  an  activation  energy  of  ss  79  kcal/mole  [13]  . A complete  review 
of  the  nitrogen- tungsten  system  is  given  by  R.  Liu  [17]  . 

* The  desorption  of  nitrogen  from  the  (110)  is  not  very  different  from 

^ that  on  other  planes.  It  is  clear,  however,  that  the  adsorption  behavior  of 

nitrogen  on  the  (110)  differs  significantly  from  that  on  the  other  planes  of 
I tungsten,  occurring  at  a much  slower  rate.  As  will  appear  in  subsequent 

chanters,  our  diffusion  studies  depend  critically  upon  this  fact. 

I 
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II.  experimental  equipment 


II. I.  Apparatus 

Prior  to  detailed  studies  of  diffusion  and  desorption  of  nitrogen  on 
the  tungsten  (110),  it  is  necessary  to  characterize  more  closely  the  inter- 
action of  nitrogen  with  the  (110)  plane  of  tungsten.  This  characterization 
is  accomplished  by  low  energy  electron  diffraction  (LEED),  flash  desorption, 
and  measurements  of  work  function  changes. 

The  experimental  apparatus,  a block  diagram  of  which  is  shown  in  Fig.  1, 
is  basically  an  old  fashioned  Varian  LEED  vacuum  system,  model  #981-0000. 

The  vacuum  system  can  be  divided  into  two  sections;  the  main  chamber,  shown 
in  Fig.  2,  in  which  all  of  the  experiments  are  performed,  and  a separate  gas 
handling  system.  A slight  modification  of  the  LEED  electronics  converts  the 
standard  LEED  system  into  a retarding  potential  analyzer  capable  of  perform- 
ing Auger  analysis.  We  also  added  a glancing  incidence  Auger  electron  gun, 
Varian  model  #981-2454.  This  gun  has  two  sets  of  electrostatic  deflection 
plates  which  make  it  possible  to  move  the  electron  beam  in  both  the  x and  y 
directions.  Figure  3 shows  a block  diagram  of  the  electronics  associated 
with  the  glancing  incidence  Auger  electron  gun.  The  Varian  Electron  Gun 
Power  module  #981-2145  provides  the  proper  potentials  and  current  to  the 
glancing  incidence  Auger  electron  gun.  The  Varian  Scanning  Samole  Positioner, 
model  #981-2157  supplies  two  sawtooth  waveforms,  which  when  app ' ied  to  the 
X and  y deflection  plates,  raster  the  electron  beam  across  the  crystal  surface. 
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Figure  1.  Block  diagram  (side  view)  of  the  experimental  apparatus 
for  diffusion  studies. 
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Figure  3.  Electronics  associated  with  the  glancing  incidence 

Auger  gun.  For  a complete  description,  refer  to  the 
text . 
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The  Scanning  Sample  Positioner  also  supplies  a signal  synchronized  to  the 
raster  waveform,  which  controls  the  raster  of  the  video  monitor.  This 
system  is  converted  into  a crude  scanning  electron  microscope  by  using 
the  current,  which  flows  through  the  crystal  to  ground,  to  modulate  the 
intensity  of  the  video  monitor. 

The  LEED  measurements  are  performed  with  a standard  Varian  4 grid  LEED 
system,  model  ?^981-0024.  Measurements  of  the  work  function  change  after 
nitrogen  is  deposited  on  the  surface  are  also  made  using  a modified  form 
of  the  same  LEED  system.  The  desorption  kinetics  of  nitrogen  from  the  (110) 
plane  are  followed  using  an  Electronics  Associates  Inc,  residual  gas 
analyzer,  model  Quad  250A. 

Connected  to  the  gas  handling  system  are  a gas  bottle,  oxygen  leak, 
Pirani  gauge,  three  nickel  getters  used  to  remove  impurities  from  the  gas 
supply,  and  two  ion  pumps.  The  manufacturer's  (Air  Reduction  Co.)  mass 
spectral  analysis  of  the  reagent  grade  nitrogen  lists  the  major  contaminants 
as:  CO  + CO2  < .5  ppm,  H < 1,  ppm,  total  hydrocarbons  < .5  ppm,  O2  < .3  ppm, 
and  argon  2.3  ppm.  The  threshold  for  other  impurities  is  < 2 ppm,  the 
balance  of  the  gas  being  nitrogen.  Further  purification  of  the  nitrogen  is 
achieved  by  selective  gettering  [4,13].  The  getter  tubes,  shown  in  Fig.  1, 
are  cooled  to  liquid  nitrogen  temperature  after  the  nickel  filaments  and 
support  loops  have  been  thoroughly  outgassed.  A current  of  2.3  to  2.5 
amperes  ac  is  passed  through  the  .010"  diameter  nickel  filaments;  at  this 
current  the  nic'xel  loops  glow  with  a red-orange  color  and  a nickel  film 
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slowly  evaporates  onto  the  getter  wall.  After  a mirror-like  nickel  coat 
develops  on  the  getter  walls  (approximately  1 hour),  the  current  is  turned 
off.  The  getters  are  then  warmed  to  room  temperature  and  all  valves  to  the 
ion  pumps  and  the  automatic  pressure  control  valve  (APC)  are  closed. 

Nitrogen  gas  is  slowly  admitted  into  the  main  section  of  the  gas  handling 

-2 

system  until  the  pressure  rises  to  5 X 10  Torr,  as  measured  by  the  Pirani 
gauge.  The  nitrogen  is  gettered  in  the  gas  handling  system  for  a minimum 
of  45  minutes  before  being  admitted  into  the  main  chamber. 

The  ultra  high  vacuum  conditions  necessary  for  reasonable  surface  studies 
are  obtained  using  a 140  liter/sec  ion  pump,  model  'r912-6006,  on  the  main 
chamber.  The  main  chamber  is  also  equipped  with  a titanium  sublimation  pump 

(TSP).  The  background  pressure  in  both  the  gas  handling  system  and  the  main 

"8 

chamber  before  bakeout  is  typically  in  the  mid  10  Torr  range.  After  a 1 
to  2 day  bakeout  at  300°C,  the  pressure  in  both  the  main  chamber  and  the  gas 
handling  system  is  generally  in  the  10  to  10  Torr  range.  A typical 
bakeout  consists  of  heating  the  system  to  300°C  for  12  to  18  hours,  cooling, 
outgassing  of  all  filaments  and  getters,  followed  by  12  to  18  hours  of 
additional  bakeout. 

II. 2.  Sample  Preparation 

The  sample  is  prepared  from  a zone  refined  single  crystal  rod,  quarter 
inch  in  diameter,  obtained  from  Material  Research  Corp.  with  the  [lOO]  along 
the  axis.  The  crystal  is  orientated  to  within  + of  the  (110)  plane  using 
Laue  back-reflection  and  then  cut  into  a wafer  .040"  thick  by  electron 
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discharge  machining.  The  sample  is  chen  mounted  onto  the  central  shaft  of 
the  polishing  jig  shown  in  Fig,  4,  and  once  again  aligned  to  + of  the 
(110)  plane,  by  Laue  back-reflection.  Bee's  wax  is  melted  around  the  rim 
and  back  face  of  the  crystal  and  the  crystal  mount  is  electrically  connected 
to  the  positive  terminal  of  a low  voltage  dc  power  supply.  The  crystal  is 
electropolished  in  a 107a  solution,  by  weight,  of  sodium  hydroxide,  using  a 
dialysis  membrane  to  separate  the  crystal  from  the  steel  cathode,  as  shown 
in  Fig.  5.  Not  shown  in  Fig,  5 is  the  tray  that  the  dialysis  membrane  and  the 
steel  plate  are  submerged  in.  The  membrane  is  approximately  one-quarter  inch 
below  the  surface  of  the  sodium  hydroxide  solution.  A potential  difference 
of  10  volts  between  the  steel  cathode  and  the  crystal  is  provided  by  a 
Hyperion  Hy-Si-20-6  power  supply.  To  polish  the  crystal,  the  polishing  jig 
is  placed  on  the  dialysis  membrane,  with  the  crystal  lightly  touching  the 
membrane,  and  moved  in  a circular  motion.  The  alignment  of  the  crystal  in 
the  jig  is  periodically  checked  by  Laue  back-reflection  and  adjusted  so  that 
the  (110)  face  is  coplanar  with  the  face  of  the  polishing  jig.  After  about 
.005"  is  removed  from  the  face,  the  crystal  is  reversed  and  the  back  face  is 
polished.  The  face  being  electropolished  is  continually  cleaned  with  acetone 
to  remove  bits  of  wax  which  become  embedded  onto  the  surface.  Both  faces 
are  repeatedly  polished  until  the  crystal  is  .015"  thick.  After  the 
polishing  is  completed,  the  central  portion  of  the  crystal  appears  flat; 
however,  it  is  very  difficult  to  prevent  preferential  etching  of  the 
crystal  edges  and  all  the  edges  show  some  rounding. 


Mounting  Rod 


Overall  view  of  the  crystal  polishing  jig.  The  tilt  screws 
are  used  to  align  the  crystal  once  it  is  on  the  mounting  rod. 
The  electrical  connection  to  the  low  voltage  power  supply  is 
on  the  back  of  a mounting  rod,  not  shown  in  this  picture. 

The  purpose  of  the  insulating  ring  is  to  electrically  isolate 
the  front  face  from  the  mounting  rod  and  prevent  it  from  being 
electropolished  while  submerged  in  the  NaOH  solution. 


Polishing  Jig 


Figure  5.  Experimental  arrangement  for  polishing  the  tungsten  crystal. 
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After  the  crystal  is  removed  from  the  polishing  jig  it  is  thoroughly 
degreased  in  acetone  and  subsequently  in  ethanol,  and  then  mounted  on  the 
holder  shown  in  Fig.  6.  This  holder  consists  of  a reservoir  attached  to 
a sapphire  rod  which  serves  as  a thermal  conduction  medium  to  the  crystal. 

By  passing  cooled  nitrogen  gas,  under  pressure,  through  the  reservoir  the 

temperature  of  the  reservoir  is  reduced  to  82°K  and  the  crystal  temperature 

O I 

to  95- K.  Connected  to  the  sapphire  rod  are  two  crystal  mounts,  made  of  i 

tungsten.  The  crystal  is  held  to  these  mounts  by  two  sets  of  tantalum 

screws  and  washers.  A tungsten  3%  rhenium/tungsten  257,  rhenium  thermocouple,  j 

.005"  in  diameter,  obtained  from  Omega  Engineering  Inc.,  is  spotwelded  to  j 

the  edge  of  the  crystal.  This  thermocouple  has  been  calibrated  against  a | 

copper-constantan  thermocouple  by  Sandstrom  and  Withrow  [l9],  so  it  may  also 

■ 

be  used  as  a low  temperature  measuring  device.  The  low  temperature  cali- 
bration curve  for  the  tungsten-rhenium  thermocouple  is  shown  in  Fig.  7. 

During  the  diffusion  and  desorption  experiments  it  is  necessary  to  i 

maintain  a high  crystal  temperature  for  extended  periods  of  time.  To 
accomplish  this  the  thermocouple  is  used  as  the  temperature  sensor  while 
the  sample  is  resistively  heated.  By  automatically  controlling  the  power 
to  the  crystal  we  are  able  to  maintain  a crystal  temperature  constant  to 
+ 3 K.  A schematic  of  the  temperature  regulating  circuit  is  shown  in 
Fig.  8.  Initially  300  amperes  ac  at  6 volts  are  needed  to  resistively  heat 
the  crystal  to  900°K.  Thereafter,  only  60-80  amperes  are  necessary  to 


maintain  this  temperature. 
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Initially  the  major  contaminant  on  the  surface,  as  shown  by  Auger 
spectroscopy,  is  carbon.  After  rasistively  hea-in?  the  sample  to  2000°K 
in  10  Torr  of  oxygen,  the  carbon  contamination  's  .educed.  Figure  9 
shows  part  of  an  Auger  spectrum  after  6 hours  of  this  heat  treatment.  After 
20  hours  of  cleaning,  the  only  detectable  contaminant  is  oxygen.  The  sur- 
face oxygen  is  easily  removed  by  heating  the  sample  to  2800°K  for  several 
minutes  once  the  ambient  gas  is  pumped  out.  An  Auger  spectrum  of  the 
clean  (110)  surface  is  shown  in  Fig.  10.  No  peaks  other  than  those  due  to 
tungsten  are  detectable. 

II. 3.  Dosing 

In  order  to  study  the  interaction  between  gas  and  metal,  a method  of 
depositing  a fixed  amount  of  gas  on  the  surface  must  be  developed.  We 
found  that  the  best  method  of  creating  reproducible  surface  coverages  is 
to  expose  the  tungsten  crystal  to  a fixed  gas  dose  from  a capillary, 

approximately  6 inches  long  and  with  a 3 mm  inner  diameter.  Dosing  is 

-2 

accomplished  by  maintaining  a constant  pressure  of  5 x 10  Torr  in  the  gas 
handling  system  and  programing  the  APC  valve  to  open  to  a predetermined 
position  for  a fixed  period  of  time  (33  seconds).  This  is  accomplished  by 
use  of  the  timing  circuit  shown  in  Fig.  il.  The  crystal  is  exposed  to 
I doses  of  two  different  sizes,  which  we  designate  as  either  large  or  small. 

I The  time  that  the  APC  remains  open  is  the  same  for  both,  but  for  a small 

i dose  the  valve  is  not  opened  completely.  An  estimate  of  the  number  of  gas 

' molecules  striking  the  surface  for  both  a large  and  small  dose  is  given  in 

I 

I Section  IV.  Although  the  absolute  value  of  the  exposure  is  uncertain  because 

of  the  possibility  of  direct  collisions  from  the  beam,  we  are  able  to  expose 
the  crystal  to  a reproducible  amount  of  nitrogen  and  obtain  reproducible 


surtace  coverages. 
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Figure  10.  Auger  spectrum  of  the  tungsten  (110)  plane  after  20  hours 
of  cleaning  in  10“^  Torr  of  oxygen  at  2000^K. 


igure  11.  Schematic  diagram  of  the  electronic  circuit  which  is 
used  to  operate  the  automatic  pressure  control  valve 
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A typical  exposure  sequence  is  initiated  with  the  cooling  of  the  sample 
to  a predetermined  temperature.  A clean  film  from  the  titanium-molybdenum 
filament  source  of  the  TSP  is  evaporated  (the  film  is  deposited  for  30 
seconds  with  a filament  current  of  45  amperes  ac),  followed  immediately  by 
heating  the  crystal  to  about  2800°K  to  remove  any  surface  contamination. 

Once  the  sample  is  thoroughly  cleaned,  as  described  earlier,  simple  flashing 
is  all  that  is  necessary  to  remove  the  surface  contamination  that  builds  up 
overnight.  The  crystal  is  then  cooled  quickly,  in  less  than  1 minute,  to  a 
predetermined  temperature  and  rotated  to  face  the  nozzle  source,  which  is 
shown  in  Fig.  2.  The  crystal  is  then  exposed  to  a dose  of  nitrogen.  During 
the  exposure,  all  filaments  are  turned  off  so  that  dissociation  of  the 
molecular  gas  by  a hot  filament  is  not  possible.  After  the  exposure  is 
complete,  the  TSP  is  once  again  turned  on  for  30  seconds  to  hasten  the  removal 
of  the  background  gas.  Taking  an  Auger  spectrum  of  the  crystal  before  and 
after  the  TSP  is  turned  on,  we  are  unable  to  detect  any  new  Auger  transitions. 
From  this  measurement  we  conclude  that  the  TSP  does  not  introduce  any 
contaminants  onto  the  crystal  surface.  During  the  exposure  period  the 
pressure  in  the  main  chamber  increased  by  about  a factor  of  5,  as  determined 
with  the  ion  gauge  during  a calibration  run,  and  not  during  an  actual 
experiment. 
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III.  FUNDA>IENTALS  A:®  TECHNIQUE 


III.l.  Auger  Spectroscopy 

To  study  the  surface  migration  of  nitrogen  we  use  scanning  Auger 

spectroscopy.  This  technique  is  ideally  suited  for  the  study  of  surfaces, 

as  it  is  extremely  sensitive,  with  the  ability  to  detect  less  chan  ,001  of 

a monolayer  '^20]  . We  define  a monolayer  as  the  number  of  atoms  on  the  top 
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layer  of  the  tungsten  (110)  plane,  which  is  1,4  x 10  atoms/cm  . Auger 
spectroscopy  is  sensitive  to  a region  extending  only  a few  atom  layers  into 
Che  crystal  [21] . The  limitation  to  the  detected  volume  is  not  the  penetra- 
tion depth  of  the  primary  beam,  but  the  escape  depth  of  the  Auger  electrons, 
which  in  general  is  10-20  A [22].  This  technique  is  therefore  sensitive  to 
changes  occurring  at  the  outer  surface. 

An  Auger  spectrum  is  most  commonly  obtained  by  the  irradiation  of  a 
material  with  a beam  of  electrons  having  an  energy  of  several  thousand  volts. 
The  Auger  process  involves  a series  of  electronic  transitions  whereby  an 
atom,  which  has  been  excited  by  an  incident  electron,  relaxes  into  its 
equilibrium  state  in  such  a way  that  it  ejects  an  electron.  The  transitions 
involved  in  such  a process  are  shown  sehcmatically  in  Fig.  12.  The  incident 
electron  creates  a hole  in  the  K level;  an  electron  from  the  L2  ^ fills 

Che  vacancy.  The  amount  of  energy  released  in  this  process  is  E - E 

^ 4,3 

This  energy  can  appear  either  as  a photon  (X-ray)  or  be  transfered  to  another 
electron.  The  latter  is  the  Auger  process  and  in  the  example  shown  in 
Fig.  12  another  L,  » electron  is  ejected  from  Che  solid.  For  any  Auger 
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Figure  12 . 


Schematic  energy  level  diagram  of  the  nitrogen  Auger 
transition. 
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transition  involving  the  energy  levels  l-fXY,  a peak  should  appear  at 


\in  ^ 


*'ana  1 


(1) 


where  E,  . is  measured  from  the  vacuum  level  of  the  analvzer,  z is  Che 

atomic  number,  A is  a correction  term  which  takes  into  account  the  change  in 

the  atomic  charge  after  Che  initial  hole  is  created,  and  5 , is  Che  work 

anal 

function  of  the  analyzer  [20],  and  not  the  work  function  of  the  crystal. 

This  can  readily  be  seen  from  Fig.  13.  Since  each  element  possesses  its  own 
"fingerprint"  of  Auger  transitions,  one  should,  in  principle,  be  able  to 
uniquely  identify  a surface  atom.  Of  course,  since  hydrogen  and  helium  each 
have  less  Chan  three  electrons  they  would  not  normally  be  expected  to  display 
Auger  transitions. 

Figure  14a  shows  the  secondary  electron  energy  distribution  as  a function 
of  emitted  electron  energy.  Near  0 volts  are  the  "true  secondaries".  Buried 
among  the  true  secondaries  and  among  secondaries  at  higher  energies  are  the 
Auger  electrons.  The  actual  number  of  these  is  small,  typically  10 
amperes  for  a primary  current  of  10  icA.  However,  since  these  transitions 
occur  at  specific  energies,  there  is  a small  but  sudden  change  in  the  energy 
distribution  curve. 

In  1968  Harris  |^21,23^  devised  a method  of  electronically  differentiating 
N(E),  the  energy  distribution  of  the  secondary  electrons  emitted  from  Che 
surface,  and  the  practicality  and  sensitivity  of  Auger  spectroscopy  for 
surface  analysis  became  obvious.  When  the  energy  distribution  curve  is 
differentiated,  the  smoothly  varying  background  remains  smooch  but  the 
differentiated  Auger  transitions  suddenly  predominate  as  illustrated  in 
Fig.  14b.  The  point  at  which  the  first  derivative  of  the  energy  distribution 
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Specimen  ap-86i 


Figure  13- 


Energy  level  diagran  of  the  complete  Auger  system 


Figure  14.  Energy  distribution  of  secondary  electrons  emitted  from 

a surface  after  electron  irradiation,  b)  the  first  derivative 
of  the  secondary  electron  distribution  in  a). 
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curve  has  its  tniniraun  is  by  custom  refered  to  as  the  Auger  energy.  The  actual 
detection  scheme  is  quite  simple,  and  is  shown  in  Fig.  15.  Our  energy 
analyzer  is  a standard  LEED  system  which  consists  of  a series  of  four 
concentric  hemispherical  grids  with  a collector.  The  first  grid  is  grounded 
to  provide  a field  free  region  between  the  sample,  which  is  located  at  the 
center  of  curvature  of  the  grids,  and  the  retarding  grids.  The  fourth  grid 
is  also  grounded  and  forms  a field  free  region  between  the  collector  and  the 
inner  two  grids. 

The  retarding  potential  analyzer  functions  by  repelling  electrons  with 
energy  less  than  E,  where  -E/e  is  the  voltage  applied  to  the  two  inner  grids. 
The  grids  act  as  a high  pass  filter  and  the  collector  current  I(S)  is 
therefore 

CS 

1(E)  a " N(S')dE'  , (2) 

“ E 

with  E'  being  a variable  of  integration.  From  Eq.  (2)  we  see  that 

I'(E)  ^ N(E)  , (3) 

where  the  ' denotes  differentiation  with  respect  to  E.  To  obtain  the 
derivative  of  the  energy  distribution  curve,  N'(E),  a small  perturbing  voltage 
AE  = k sin(a)t)  is  applied  to  the  two  inner  grids.  The  collector  current  1(E) 
is  thus  energy  modulated.  Expanding  I(E  4-  AE)  in  a Taylor  series,  we  find 

2 3 

I(E  + AE)  = 1(E)  + I'(S)iE  + + ...  . (4) 
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upon  subscitution  and  rearrangement  this  becomes 

3 

T"  'k 


I = Iq  + .r'k  + -3 


2 4 

T'’V 

+ . . .]sin('jut)  - [— ^ + + . ...cosfxt)  , (5) 


or 

2 

r = I(Eo)  + I'k  sin(xt)  - cos(2xt)  , (6) 

3 

provided  k and  higher  terms  can  be  neglected. 

3y  using  phase  sensitive  detection,  M(E)  can  be  obtained  on  tuning  to 
the  modulating  frequency;  detection  at  2;ju  gives  N'(E),  which  is  commonly 
refered  to  as  the  Auger  spectrum.  A complete  review  of  Auger  spectroscopy 
is  given  by  C.  C.  Chang  [20], 

To  electronically  differentiate  the  energy  distribution  we  modified  a 
standard  4 grid  LSED  system,  as  shown  in  Fig.  15,  for  use  as  an  energy 
analyzer.  The  purpose  of  the  capacitance  neutralization  circuit  (CNC)  and 
the  active  filter,  shown  in  Figs.  16  and  17,  are  to  attenuate  that  part  of 
the  collector  current  whose  frequency  component  is  the  same  as  the  modulation 
frequency  x and  to  permit  only  the  second  harmonic  of  the  modulation  frequency 
(2x  = 4100  Hz)  to  enter  the  lock-in.  The  CMC  forms  a capacitance  bridge 
between  the  inner  two  grids  and  the  collector  on  one  hand,  and  the  collector 
and  ground  on  the  other.  By  properly  matching  the  capacitance  of  the  CNC 
to  the  LEED  system,  the  coupling  to  the  collector  is  reduced,  minimizing  the 
collector  current  induced  at  the  modulating  frequency  x.  This  prevents  the 
Ithaco  391  -A  lock-in,  used  in  our  experiments,  from  overloading  on  the 
fundamental  of  the  perturbation  frequency  and  permits  us  to  examine  that 
part  of  the  collector  current  that  goes  as  2x.  The  active  filter  not  only 
attenuates  the  fundamental  but  also  reduces  the  bandwidth  and  increases  the 
gain  of  the  system  to  give  an  order  of  magnitude  improvement  in  the  signal- 


to-noise  ratio. 
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Figure  16.  Schematic  of  the  capacitance  neutralization  circuit 
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Of  the  three  najor  sources  of  noise  [24],  flicker  noise,  thermal  noise, 
and  shot  noise,  shot  noise  is  the  largest  that  cannot  be  reduced.  The  shot 
noise,  I , in  thermionic  circuits  is  given  by  ’25] 


= 2el  R 
n 0 


(7) 


where  e is  the  electronic  charge,  the  total  collector  current  and  3 the 

bandwidth.  Under  typical  operating  conditions,  with  = 100  -A,  3 = . 04 

(for  a time  constant  of  3 sec  and  a 12  d3/octave  rolloff),  the  rms  noise 

voltage  is  =3  1 uV,  From  a clean  surface  the  signal  voltage  S can  be  larger 

than  1000  uV  _22-24j  so  'that  the  signal-to-noise  ratio  is  S/N  = 1000.  The 

Auger  signal  can  originate  several  atom  layers  beneath  the  surface,  this 

15  2 

means  that  the  signal  is  caused  by  as  5 x 10  a toms /cm  . If  the  limit  of 

12  2 

detection  is  defined  as  S/N  = 1,  then  as  little  as  5 x 10  atoms/cm  can 

be  detected.  The  only  way  to  increase  the  sensitivity  is  to  increase  or 

to  decrease  3.  Since  the  signal  is  proportional  to  and  the  noise  is 

proportional  to  v^I^,  hy  increasing  the  3/N  ratio  increases.  The  S/N 

ratio  can  also  be  increased  by  increasing  the  collection  efficiency. 

Decreasing  3 will  also  improve  the  S/N  by  decreasing  the  noise  current  I ; 

n 

however,  this  is  unattractive  since  it  will  require  an  increase  in  the 
already  long  time  constant. 

In  some  cases  .Auger  spectroscopy  can  be  used  to  obtain  quantitive 
Information  as  to  the  surface  concentration  of  a particular  element, 
relative  to  3 standard  concentration.  That  this  is  the  case  can  be  seen 
by  examining  the  energy  distribution,  N/F),  for  an  Auger  transition.  The 
energy  distribution  of  most  .Auger  transitions  can  be  described  by  a 
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Gaussian  or  Lorentzian  curve.  If 

, for 

example,  we 

take  the  energy  dis- 

( 

». 

'' 

1' 

tribution  to  be  a Gaussian,  as  in 

F ig. 

13a,  then  I' 

(E)  is  of  the  form 

► 

i 

I'(E)  = i;- 

1 

' i 

.,2  _ 

exp  1 - 1 

, (8) 

' 

/2m  ■' 

- 2c~  - 

1 vhere  2c  is  the  full  width  at  half  height  of  the  curve  in  Fig.  13a.  The 

peak  anplituda  of  the  Auger  spectrun  is  given  by  *24] 


A^  = ± .06  .|/i  , (9) 

where  k is  the  amplitude  of  the  perturbation,  i is  the  total  Auger  current  of 
a single  transition  and  z corresponds  to  the  inflection  point  of  the  N(E) 
cujr/e.  In  this  case  the  Auger  peak-to-peak  height  is  proportional  to  the 
.Auger  current  i.  Two  points  must  be  kept  in  mind  before  the  differentiated 
peak-to-peak  height  can  be  used  as  a measure  of  the  surface  concentration; 

1.)  The  differentiated  peak  must  be  symmetric,  as  it  is  for  a Gaussian, 
Lorentzian  or  other  reasonable  Auger  electron  energy  distribution;  this  can 
be  seen  in  Fig.  ISb.  2.)  The  A-iger  peak  width  must  remain  constant, 
independent  of  the  surface  coverage.  An  increase  in  the  surface  concentration 
will  increase  i and  not  c.  Using  Eq.  (8)  we  see  that  the  peaks  on  the  Auger 
spectrum  will  always  appear  at  ±c,  independent  of  the  value  of  i.  If 
these  two  conditions  are  not  met  than  the  peak-to-peak  height  of  an  Auger 
transition  will  not  be  indicative  of  the  relative  surface  coverage. 

Fortunately  these  conditions  are  approximately  satisfied  by  most  elements. 
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III. 2.  Surface  Characterization 

III. 2. A.  Surface  Scructure 

There  is  no  technique  available  for  accurately  determining  the 
perfection  of  a surface.  Gross  atonic  disarrangements  can  be  detected 
by  electron  microscopy;  however,  surface  defects  smaller  than  ss  500  A 
cannot  easily  be  detected.  LEED , on  the  other  hand,  can  provide  informa- 
tion only  about  surface  structures  which  are  periodic  within  a region  over 
which  the  primary  wave  is  coherent.  The  coherence  length  is  limited  by; 
(a)  spatial  coherence,  which  is  dictated  by  the  angular  divergence  of  the 
primary  beam,  and  (b)  temporal  incoherence  which  results  from  the  energy 
spread  of  the  primary  beam.  The  coherence  length  is  related  to  the 
angular  divergence  3^  and  the  energy  spread  AE  of  the  primary  beam  by 
[26] 


L 

s 


2(1#33 


(10) 


where  X is  the  electron  wavelength  and  S is 
This  gives  a coherence  length,  L,,  of  about 
beam  energy  of  100  volts. 

In  order  to  understand  the  nature  of  a 
to  characterize  the  degree  of  perfection  of 


the  normal  electron  energy. 
2 3 » 

10  to  10  A,  tor  a primary 


surface  layer  it  is  necessary 
the  surface  as  well  as  the 


structure  of  the  adsorbed  layer.  The  periodicity  of  an  adsorbed  layer 
can  be  characterized  by  the  use  of  low  energy  electron  diffraction,  LEED. 
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To  derive  the  iaterferer.ee  function  ve  describe  the  scattered  amplitude 
for  a 2 -d itnens iona  1 array  of  scatters  by 


n 


f 

n,!n 


e.xp  (!.!:<' r ) 
n 


(11) 


where  2.K  = being  the  scattered  wave  vector,  the  incident 

wave  vector,  f the  atomic  scattering  factor,  and  r the  position 
n,m  = n,m  ^ 

vector  of  t.he  atom  located  at  nax  + mby  (n  and  m being  integers),  wit.h  a 
and  b the  lattice  constants  as  shown  in  Fig.  19.  The  scattered  wave 
vector  K,  can  be  written  as  [27] 


-<-  = (cos  a)x  + — (cos  3)y  + n—  (cos  y)2, 

• K K \ ' 


In 


(12) 


where  a,  3 and  V are  the  angles  formed  by  and  the  x,  y and  2 axis 
respectively.  For  normal  incidence. 


AK  = K^*r 

t n,m 


(13) 


Equation  (11)  then  becomes 


■J  * V-  Z f expr  - — — ^ (r.acosc2  - mbcosS  )'^ 
n.m  ^ ' 


The  interference  function,  F,  is  the  square  of  'It  or 


2**  -I ' 2 

F « ^ Z f cos[-r^  (nacoscy  -i-  mbeosS)!,* 
■a,m  ^ 

' _ - - 2-t  . - ' 2 

+ 2,  t sini.-m-  (nacosa  -r  mocosc)J'-  , 

V n , m A. 

a,m 


(14) 


(15) 


Figure  19.  ?(2;<2)  domains  on  the  (ilO)  plane  of  tungsten;  d, 

the  translation  vector  between  the  two  domains.  ^ 


J 
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where  coso;  = (— )h,  c033  (r)k,  '/■.  Che  elecCron  wavelength  and  h,k  che  indices 

3 0 

of  the  diffracted  beam.  By  calculating  the  value  of  F for  a particular 
surface  structure  we  obtain  Che  syrnciecry,  spacing  and  beam  profiles  of  Che 
diffraction  pattern,  but  not  their  intensities. 

Under  certain  circumstances,  however,  additional  information  about  the 
atomic  arrangement  can  be  obtained.  An  adsorbed  gas  sometimes  forms  domains. 
Vhen  the  dom.ains  are  out  of  phase  with  one  another,  splitting  of  some  of  the 
beams  may  occur.  This  can  be  seen  by  examining  the  scattering  from  two 
domains,  separated  by  a distance  (dj^),  as  shown  in  Fig.  19.  Using  Sqs.  (11) 
and  (13),  the  scattered  a.mplitude  for  che  first  domain  is  given  by 


- 2 t expfiKgT  j 

1 n,m  ^ f n,m 

n ,m 


(16a) 


and  for  the  second  domain  denoted  by  2, 


'9  ^ axp[;i:<--r  , ,]  . 

^ —I— rn,[n 

n ,01 


(15b) 


, 2 

The  interference  function  for  both  dom.ains  (~^2^  given  by  ['!f,  , 


1 


that  is  by 


F„  = c.  f ex?(iX.-r  ) + T f ,exp(iK,-r  , ,)] 

12  ^ n,m  t n,m  , , n ,m  f n ,m 

n,m  n' ,m' 


(17a) 


F,o=FiTF_+i-  ^ r r,  , [exp(r.<.T  )exo(-i.<^-r  , ,) 

12  1 2 , , n,m  n',m'  ^ t n,m  f n ,m 

n,m  n ,m 


X exp(-iK,T  )exp(iK-'r  , ,)]  , 


(17b) 


where  F^^  = and  F2  = If  the  domains  are  of  equal  size,  Chen 
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Fi2  = 2?^  +2  2 cos  (K.-d^)  , (17c) 

‘ n ,ni 

where  is  the  translation  vector  between  the  tw’o  domains.  Equation 
(17c)  can  be  simplified  to  [25] 

Fi2  = 2F^[1 +cos(K--d2)]  . (17d) 

Kj 

The  maxima  of  F occur  when  = ha'"  + kb"  , where  k and  h are  the  indices 
of  the  diffracted  beam  and  a"  and  b*  are  the  reciprocal  lattice  vectors, 
vvhen  one  of  these  maxima  coincides  with  a minimum  of  the  (1  + cos  (K - -d  j^) ) 
term,  beam  splitting  occurs.  This  happens  when 

hn  + km  = y (2i  -t  1)  , i = integer  , (18) 

where  dj^  = na  -r  mb. 

3y  application  of  Eq.  (15)  to  a proposed  surface  structure,  the  symmetry 
spacing  and  beam  profile  of  the  resulting  LEED  pattern  can  be  calculated. 

Beam  splitting  m.ay  be  used  to  obtain  insight  into  the  growth  mechanism  of 
the  adlayer,  its  size,  the  direction  of  the  boundary  and  how  well  formed 
the  domains  are. 

111.2.3.  Surface  Defects 

Perhaps  the  most  common  surface  defects  are  steps  or  terraces.  Ellis 
and  Schwoebel  [29]  and  Kenzler  [30]  have  shown  theoretically  and  experi- 
mentally that  a regular  array  of  surface  steps  causes  a characteristic 
splitting  of  many  of  the  LEED  beams  into  doublets  or  even  multiplets.  The 
terrace  width  can  be  determined  from  the  separation  of  the  spots  which  make 
up  the  doublet. 
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A diagram  of  a surface  vith  step  edges  running  along  the  x axis  is 
shotTO  in  Fig-  20;  Fig-  21  shows  the  projection  or  this  plane  along  the 
step  direction.  If  we  consider  scattering  from  only  the  top  most  layer 
of  Figs.  20  and  21  and  use  Eq.  (11),  the  scattering  factor  for  the  first 
row  of  N -f- 1 lattice  sites  along  the  y direction  is 

N 

= 'if  (0,0)  = f I ex?[iqK.-a]  , (19) 

y y q=0 


where  f is  the  atomic  scattering  factor  of  each  surface  atom.  For  the 
second  row  the  scattering  factor  is  given  by. 


N _ _ ^ 

jf  (b,0)  = f Z eXD [iK- • (aa +b)  ] 
y q=0  ‘ " 


(20) 


The  scattering  factor  from  the  first  terrace  of  M-i-l  rows,  where  M is  odd,  is 

r ^ ^ N ^ ^ ^ N _ ^ ^ 

'*  = f i 2)  exp(iK."qa)  + I exp[iKp' (qa +b)  ] + H exp[iK,' (qa +2b)  ] + . - - j- 

q=0  “ q=0  q=0 


N 


t(M-l) 


t [ 1 + exp (l.<- .b)  ] i.  exp[].!C,-qa]  E expriKc-2sb]  . (21) 

^ n t n ‘■ 

q=0  3=0 


The  scattering  factor  from  subsequent  terraces  is 


.j,(2)  _ .1,(1)  exp[iK£‘(Na  + s^^  + S2  + s^)] 


exp[iX^-(2Na  + 2^^  ^27^)]  . 


(22) 


(T) 

The  total  scattering  factor  if  of  n steps  and  n-i-l  (n  odd)  terraces  is 


,^(T)  . .,/p) 

p=l 


AP- 3^J 


Figure  21 


I 


Projection  of  the  stepped  surface  of  Fig.  20  along  the 
step  direction. 
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= [ 1 + exp (iK . • (Na  + s,)]  X exp[iK.-(Ma  - + 3^)2t] 


Equation  (23)  can  now  be  '.•rritten  as 


rr)  ^ - t(M-l) 

= f [ 1 -f  exp(iK- •a')  ] L exofiK^-aa]  I exp[iK-'2sb] 
^ ' q=0  ■ ^ ‘ s=0  ■ 


X [ 1 T exp  (iK  - ■ (2Ia  -’r 


t(M-l) 


2 'y  -rcr 

=1  =2 


E exp[i:<..(Na  + g g^)2t] 
:=0 


The  intensity  distribution  is  given  by  the  square  of  the  absolute  value 
of  the  scattering  factor, 


■>  2 — — ' 2 r.'’  ""  ”*  ''■"'i' 

I = 4f"  cos  (K,-a)*i4  cos  i ■ vf- ; • v Na  +g,  Tg„  +g.  ■'  ' 
I ' _ ' z \ i l }>  ■ ■ 


sin2^r^^i,.k 


-K.-a- 


!in^;  sin^M'SlljK,-  (Na  ■ 


sin" [K. -b] 


sin"[K^-  (Na  +g^  )] 


Transforming  this  into  the  coordinate  system  illustrated  in  Fig.  22,  where 


sinr  sin  X 


"Z-bl 

■ — ; — / cos-! 

\ / 


, 2*^  \ 

— : — . sin9  cos  x 


- - /^"Sp' 

Kf-g,  = sin9  cos  X 


:<f -33 
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4P-9U 


Figure  22.  Geometry  appropriate  to  scattering  from  stepped  surface. 

Terraces  are  parallel  to  the  xy  plane  and  steps  are 
parallel  to  the  x axis. 
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Therefore,  the  intensity  at  an  angle  X,  when  9 = 90  is  given  by 


I(X)  = const. 


sin“i  sinX  ' sin^[(n+l)  j (Na  sinX  +g^  sinX  + (1  + cos  X ) ) i 


. zrxa  . ^ 1 
sin  ' -IT  sinX  I 

Li  -J 


sin*^  [K(Na  sinX  +g  sinX  +g,  (1  +cos  X ))] 


where  K = ^ • The  fir^t  term  corresponds  to  diffraction  from  the  terrace. 
The  main  maxima  occur  at  angles  X,  where  ^ Ka[sin(X)]  = nn.  In  the  limit  of 
an  infinite  number  of  steps,  and  using  the  relationship 


2 

sin  [ (n+1)  x]  => 


= I 6(x-iTT)  , 


n— 00  s in  X 


the  second  term  can  be  replaced  by 


* 1 
2 0[j  K (Na)sin  X + j fCg^  sin  X + 2 Kg2  ^ 

i=-® 

This  term  depends  only  on  the  step  width  Na , the  step  height  g2 , and  hori- 
zontal displacement  gj^.  The  separation  of  two  adjacent  delta  functions  is 


now  given  by 


[(Na)cos  X- Sj^sinX  +g  cosx  ] ’ 


and  for  X 0, 


Na  +g2 

From  this  expression  the  terrace  width  can  easily  be  determined  by 
measuring  the  angular  separation  of  the  (0,0)  doublet  at  near  normal  incidence. 
From  a practical  point  of  view,  this  technique  can  only  be  used  if  the  step 
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■ 9 

I width  is  less  than  about  150  A-  For  steps  wider  than  this  the  beam 

I Splitting  is  so  small  that  our  LEED  system  is  incapable  of  resolving  the 

doublet. 

H 

III.2.C.  Measurements  of  Work  Function  Changes 

The  difference  in  the  electrostatic  potential  between  the  inside  and 
outside  of  a metal  is  termed  the  surface  potential.  The  related  properties- 
the  work  function  and  the  contact  potential  difference,  measure  the  free 
energy  change  that  occurs  when  electrons  are  moved  from  one  conductor  to  a 
vacuum  and  from  one  conductor  to  another,  respectively.  The  formation  of 
an  adlayer  modifies  the  dipole  structure  of  the  metal  surface  and  therefore 
changes  the  work  function.  The  effect  may  be  evaluated  in  terms  of  the 
dipole  moment  produced  by  the  adlayer.  It  can  be  shown  that  the  change  in 
the  value  of  the  work  function  corresponding  to  the  formation  of  a layer 
with  a fractional  coverage  9 is 

' A(2)  = 4r:s  MG  (32) 

: ^ 

where  M is  the  dipole  moment/atom  and  s the  total  sites  available  per 
I m 

unit  area  of  the  adsorbate.  A review  of  techniques  used  for  measuring 

surface  potentials  is  given  by  Culver  and  Tompkins  [31]  and  by  Knapp  [32]. 

, To  study  the  change  in  work  function  after  the  adsorption  of  a gas, 

our  LEED  system  is  modified  to  be  used  as  a retarding  diode  analyzer. 

\ 

The  retarding  diode  technique  is  a variation  of  the  contact  potential 
' method  described  by  Langmuir  and  Kingdon  [33]  and  later  by  Skelton  [34], 

The  only  equipment  needed,  in  addition  to  the  standard  LEED  system,  is  a 
microammeter  and  a small  dc  power  supply.  The  experimental  arrangement 
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Figure  23.  Block  diagram  of  the  retarding  diode  analyzer  for  measuring 
work  function  changes. 


Current  (xlO"^Amp) 
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I 

Figure  24.  Characteristic  dc  calibration  curves  obtained  at  different 
retarding  voltages.  In  this  case  a 1 volt  change  in  the 
bucking  voltage  produces  a 0.95  volt  shift  in  the  I-V  curve, 
a)  bucking  voltage  = 20  volts,  b)  bucking  voltage  = 20.5 
volts,  c)  bucking  voltage  = 21.5  volts;  actual  shift  in 
I-V  curve  0.95  volt. 
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is  shown  in  Fig.  23.  Retarding  potential  curves  are  derived  by  measuring 

the  current,  I , to  the  crystal  as  the  accelerating  voltage  V is  increased, 
a 

In  our  system  the  primary  beam  voltage  from  the  LEED  gun  is  set  at  21 

-7  -8 

volts  and  the  current  to  the  crystal  is  between  10  to  10  amperes.  To 
check  the  accuracy  of  the  retarding  potential  analyzer,  a dc  calibration 
is  performed.  Initially  the  primary  emission  current  is  ^ 50  uA  and  the 
accelerating  voltage  and  the  bucking  voltage  are  identical.  The  accelerating 
voltage  is  then  slowly  increased.  At  some  low  voltage  (relative  to  the 
bucking  voltage)  the  electrons  have  barely  enough  energy  to  arrive  at  the 
crystal.  As  the  acceleration  voltage  V is  increased,  the  current  to  the 
crystal  increases  and  eventually  saturates  (typically  at  about  5 x 10  ^ 
amperes).  This  is  then  repeated  for  a number  of  different  bucking  voltages, 
and  a dc  calibration  is  obtained.  A sample  calibration  is  shown  in  Fig.  24. 
There,  a 1 volt  change  in  the  bucking  voltage  produces  a 0.95  volt  shift 
in  the  I-V  curve;  the  gain  of  the  system  is  therefore  0.95.  The  shift  of 
each  I-V  curve  to  the  right  with  increasing  bucking  voltage  is  due  to  the 
change  in  potential  of  the  crystal;  if  the  surface  potential  itself  changes, 
the  I-V  curve  would  also  shift,  and  this  shift  would  be  proportional  to 
the  change  in  work  function.  In  our  experiment  the  accelerating  voltage  is 
held  constant  and  instead  the  bucking  voltage  is  varied;  the  accelerating 
voltage  is  always  kept  at  21  volts.  The  difference  in  voltage  between  the 
bucking  voltage  and  the  accelerating  voltage  is  the  retarding  voltage  in 
the  I-V  curve. 

Consider,  for  example  what  happens  when  the  work  function  increases 
after  a gas  is  deposited  on  the  surface.  At  any  particular  bucking  voltage 
V,  the  current  decreases  and  the  I-V  curve  shifts  to  the  right.  If  the  slope 
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of  the  I-V  curve  (in  the  linear  region)  is  constant,  the  change  in  I is  proper 
tional  to  the  change  in  work  function  and  the  proportionality  constant  is  the 
slope  of  the  I-V  curve. 

When  the  LEED  system  is  used  in  this  manner  to  measure  work  function 
changes,  there  are  three  possible  sources  of  error.  The  first  is  a change  in 
the  work  function  of  the  filament.  Since  the  gun  is  not  operated  in  the 
space  charge  limited  mode,  the  adsorption  of  nitrogen  on  the  cathode  may 
affect  the  changes  observed  in  the  work  function  of  the  sample.  Fortunately, 
the  temperature  of  the  filament  is  high  enough  that  any  adsorbed  nitrogen 
will  be  desorbed  at  the  operating  emission  current.  In  addition  it  can  be 
shown  that  the  emission  current  to  the  anode,  I^,  is  given  by 

where  T.  is  the  filament  temperature,  A a constant,  the  work  function 

of  the  anode  and  is  the  difference  in  potentials  between  the  cathode  and 

the  anode.  From  Fig.  25  we  see  that  the  anode  current  in  the  diode  is  not 

dependent  on  the  value  of  the  cathode  work  function,  provided  that 

< [ 'i  . -V  ],  where  5 is  the  work  function  of  the  cathode, 
anode  a cat 

The  second  possible  source  of  error  lies  in  a change  in  the  reflection 
coefficient  r when  gas  is  adsorbed  on  the  surface.  Equation  (33)  properly 
should  include  a factor  (1-r),  where  r is  the  reflection  coefficient  of  the 
electron  at  the  anode  and  is  taken  to  be  independent  of  electron  energy. 

A change  in  the  reflection  coefficient  would  manifest  itself  as  a change  in 
the  slope  of  the  I-V  cur’/e  with  surface  coverages.  Defocussing  of  the 
electron  beam  constitutes  yet  a third  source  of  error.  For  small  values  of 
I^,  the  electrons  arriving  at  the  crystal  have  almost  zero  energy. 
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Figure  25.  Energy  level  diagram  of  the  retarding  potential  diode. 
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The  electron  beam  is  thus  easily  defocussed  and  more  of  the  electrons  miss 
the  sample.  Beam  defocussing  is  reduced  by  using  a large  bucking  voltage 
and  by  operating  the  gun  in  such  a way  that  electrons  leave  the  gun  at 
high  velocities  and  decelerate  only  when  they  arrive  at  the  crystal.  This 
is  accomplished  by  biasing  the  drift  tube  of  the  LEED  gun  as  shown  in 
Fig.  23. 

The  usefulness  of  surface  potential  measurements,  using  the  retarding 
diode  technique,  lies  in  its  simplicity,  minimal  eq  iment  requirement  and 
sensitivity  to  changes  that  occur  on  the  surface.  The  retarding  diode 
technique  can  easily  detect  changes  of  less  than  0.01  volt  in  the  surface 
potential.  Since  the  adsorption  of  a monolayer  on  the  surface  typically 
produces  a change  of  several  tenths  to  one  volt  in  the  surface  potential, 
small  changes  in  surface  coverage  should  readily  be  detected. 
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IV.  RESULTS 

IV. 1.  Specimen  Characterization 

IV.  1. A.  Contamination 

The  major  sources  of  surface  contamination  are  impurity  diffusion  from 
the  bulk  and  sample  holder,  from  the  background  gas  and  possibly  from  the 
nitrogen  gas  supply.  To  insure  the  purity  of  the  nitrogen  prior  to 
exposure  to  the  crystal,  the  gas  is  passed  over  freshly  evaporated  nickel 
films  before  the  start  of  an  experiment. 

Prior  to  an  experiment  the  sample  is  flashed  to  2800°X  and  an  Auger 
spectrum  of  the  tungsten  surface  is  taken  to  insure  its  cleanliness. 

The  crystal  is  then  exposed  to  nitrogen  and  another  .‘^uger  spectrum  is  taken. 
The  only  new  peaks  visible  are  the  nitrogen  triplet  which  occurs  at  about 
380  volts;  this  is  shown  in  Fig.  26,  From  this  we  conclude  that  any 
contamination  from  the  gas  supply  is  too  small  to  be  detected.  The  only 
contamination  found  results  from  the  dissociation  of  the  background  gas 
on  the  surface,  kfter  4 to  5 hours  in  a vacuum  of  10  to  10  Torr,  as 
measured  by  the  nude  ion  gauge  in  the  main  chamber,  a slight  carbon  peak 
is  found,  amounting  to  approximately  one-half  that  of  the  carbon  peak  in 
Fig.  9.  The  surface  carbon,  presumably  arising  from  the  adsorbtion  of  CO, 
is  easily  removed  by  a flashing  to  several  thousand  °K  for  a few  seconds. 

IV.l.B.  Surface  Arraneement 

Perhaps  the  most  common  surface  defects  are  steps.  In  Sec.  III.2.B. 
we  discussed  a method  by  which  LEED  measurements  are  used  to  determine 
the  step  density.  At  several  different  locations  along  the  crystal  edge 
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Energy  (eV)  ^ 

Figure  26.  Auger  scan  of  the  tungsten  (110)  plane  after  nitrogen 
deposition  on  the  surface.  The  peaks  due  to  nitrogen 
are  located  at  » 380  volts. 
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we  observed  evidence  of  beam  splitting  due  to  surface  steps.  This  is  evident 
in  the  LEED  picture  of  Fig.  27,  which  shows  the  (0,0)  beam,  slightly 
displaced  from  the  center  of  the  LEED  screen,  in  order  to  be  photographed. 

The  primary  beam  energy  is  35  volts  and  from  the  splitting  we  are  able  to 
estimate,  using  Eq . (31),  that  near  the  crystal  edge  the  terrace  widths 
are  approximately  30  atoms  wide.  The  step  edges  are  found  to  be  parallel 
to  the  [113]  direction.  When  the  LEED  beam  is  moved  to  the  center  of  the 
crystal  no  splitting  is  evident.  This  is  not  meant  to  imply  that  there 
are  no  steps  in  the  center  of  the  crystal.  In  the  center  the  steps  could 
be  arranged  in  a random  array,  or  the  terraces  may  be  wider  than  about 
30-40  atoms  spacings.  For  terrace  widths  exceeding  30-40  atom  spacings 
the  angular  separation  of  the  doublet  is  less  than  1°  and  we  would  be  unable 
to  resolve  such  a doublet. 

i 
1 

IV. 2.  Ordinary  Nitrogen  Lavers  on  W(llO) 

IV. 2 -A.  Adsorption  and  Desorption  of  v Nitrogen  from  the  W(llO) 

Since  we  are  using  electron  stimulated  dissociated  of  y nitrogen  to 
create  our  localized  surface  deposits  (as  will  be  discussed  in  greater 
detail  in  Sec.  IV. 3. A),  we  need  to  better  understand  the  interaction  of  V 
nitrogen  with  the  (110)  plane  of  tungsten.  To  this  end  flash  desorption 
is  used  to  follow  the  desorption  kinetics  of  y nitrogen.  After  the  crystal 
is  cooled  to  95°K  and  exposed  to  predetermined  amounts  of  nitrogen,  the 
crystal  is  resistively  heated  to  a temperature  T « 300°K.  A residual  gas 
analyzer  monitors  the  amount  of  nitrogen  liberated  during  the  heating 
interval,  as  a function  of  T.  Both  mass  14  and  28  peaks  are  compared  and 


Figure  27. 


Splitting  of  the  (0,0)  beam,  due  to  a regular  array  of  surface 
steps,  near  the  edge  of  the  crystal.  The  primary  beam  voltage 
is  35  volts. 
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the  desorption  rate  curves  are  identical  in  shape  and  peak  desorption 
temperature.  It  is  concluded  from  this,  and  the  surface  contamination 
measurements  mentioned  earlier,  that  the  contamination  from  CO  (also  mass  28) 
is  insignificant. 

In  general,  the  rate  of  thermal  desorption  of  a gas  from  a unit  surface 
area  is  given  by  [35,36] 

N(t)  = - = V C''(t)exp[-E  /RT]  , (34) 

at  n aes 

where  N is  the  number  of  atoms  desorbed  per  unit  area  per  unit  time,  n the 
desorption  order,  the  characteristic  frequency,  the  activation  energy 

of  desorption,  and  C(t)  the  surface  concentration.  In  flash  desorption 
measurements,  a surface  containing  an  adsorbate  at  an  initial  concentration 
Cq  is  rapidly  heated,  usually  in  such  a way  that 

T = Tq  + a t , (35) 

Tq  being  the  initial  temperature  and  oc  the  heating  rate.  The  pressure  rise 
in  the  system  is  proportional  to  if  the  pumping  speed  is  hign.  This 

pressure  rise  can  be  recorded  as  a function  of  temperature  or  time.  By 
integrating  the  desorption  rate,  Eq.  (34),  the  number  of  surface  atoms 
desorbed  can  be  calculated.  This  is  accomplished  by  measuring  the  area 
under  the  flash  desorption  curve.  By  exposing  the  crystal  to  increasing 
doses  of  nitrogen,  followed  by  the  integration  of  the  obser/ed  flash  desorp- 
tion curve,  the  surface  coverage  per  dose  can  be  calculated,  relative  to 
some  standard  exposure , For  small  doses  the  coverage  is  found  to  be  propor- 
tional to  the  exposure.  Figure  28  shows  the  amount  of  nitrogen  desorbed  for 
exposures  ranging  from  one  to  six  small  doses.  All  areas  under  the 
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desorption  curve  are  relative,  with  the  area  for  a saturated  surface  taken  as 
unity.  For  large  doses,  the  amount  of  nitrogen  desorbed  became  non-linear  for 
coverages>  0.4  of  the  saturated  surface.  Figure  29  shows  the  relative  amount 
of  nitrogen  desorbed  after  exposures  of  1 to  6 large  doses.  If  we  assume 
that  the  saturation  coverage  of  Y nitrogen  on  the  (110)  plane  is  one  mono- 
layer  or  1.4  X 10^^  molecules/cm^ , and  the  sticking  coefficient  for  small 

coverages  is  ~ 1 [37],  then  one  small  dose  corresponds  to  an  impingement  of 

13  2 

about  7 X 10  molecules/cm  , and  one  large  dose  is  equal  to  an  i.mpingement  of 

14 

of  5.88  X 10  molecules /cm” . 

Figure  30  shows  several  thermal  desorption  curves  plotted  on  the  same 

o 

scale.  Within  the  large  peak  centered  near  152  K,  there  appear  to  be  at  least 

1,  and  probably  more,  discrete  sub-peaks.  The  major  peak  is  associated  with 

V nitrogen  desorbing  from  the  (110)  plane,  the  smaller  one  (ones)  at  higher 

temperature  are  possibly  due  to  desorption  from  the  edges.  The  crystal  is  a 
1 " 

disc  ~ ^ in  diameter  and  0.015"  thick.  This  results  in  an  edge  area  which 
is  117o  of  the  total  crystal  area.  Therefore,  some  of  the  desorbed  nitrogen 
originated  from  planes  other  than  the  (110) . At  still  higher  temperatures 
(above  1400°K)  we  notice  the  desorption  of  3 nitrogen. 

For  a first-order  desorption  process  the  temperature  T^  at  the  peak  in 
the  desorption  curve  can  be  found  from  Eq . (34),  to  be 

%.4e,p(-E  /RT  ) . (36) 

RT 

P 

The  shape  of  the  entire  desorption  curve  can  also  be  calculated  using  Eq.  (34), 
assuming  a high  pumping  speed.  For  a linear  heating  schedule  and  a first 
order  desorption  process,  integration  of  Eq.  (34)  yields 


Desorption  Rote  molecules  (cm’’^sec 


63 
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Temperature  (®K) 


AP-830 

Figure  30.  Desorption  rate  of  molecular  nitrogen  from  the  (110)  plane 
of  tungsten.  The  nitrogen  is  initially  deposited  at  95°K. 

The  three  initial  coverages  are  1.0,  0.32  and  0.45  monolayer. 
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r dC 

- c 


1.1 

1 -> 


= T j exp(-E  /RT)dT  . 


= In  ' for  n=l, 


where  is  the  coverage  at  the  start  of  the  temperature  sweep,  and  C is 

u p 

the  coverage  at  T^.  Substituting  u = -E^^^/RT,  integrating  by  parts  and 
using  the  identity 


exp(t)dt  exp(u)  r,  , fiW  . fzW  . 

-CO  y 


for  u » 1 , (38) 


we  obtain  the  general  shape  of  the  flash  desorption  curve  as 


N _ I :“^des  fl  l"\  , fT 
p p p 


^des  ;l  _ 


R M T 

P 


C39a) 


-2  -1 

Here  N is  the  number  of  desorbed  species  (cm  sec  ),  and  N the  number 

P 

-2  -1 

desorbed  (cm  sec  ) at  T^.  For  a surface  species  bound  in  j different 
states,  the  shape  of  the  desorption  curve  for  (n=l)  is  given  by 


Total 


des  f _1 
R Vt 


"des  f\  1 


P 


4 - ) -1  V 

^(i)  h 


(39b) 


Even  without  a complete  analysis,  a number  of  important  facts  can  be 
derived  for  the  desorption  kinetics  of  nitrogen  from  the  (110)  plane  of 

V 

tungsten  [36,39].  The  peak  desorption  temperature  of  the  v states,  T^,  is 
the  same  independent  of  initial  coverage.  This  can  only  be  true  for  a 
first-order  desorption  process  with  a constant  desorption  energy. 
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Looking  at  the  high  temperature  side  of  the  y desorption  curves  we  see  that 


the  curves  have  a shoulder.  This  is  probably  due  to  another  molecular 

entity  desorbing  with  a slightly  different  energy.  If  we  assume  that  there 

is  only  one  state  on  the  high  temperature  side,  and  that  it  follows 

first-order  desorption  kinetics,  then  we  can  easily  separate  the  two 

desorption  processes.  This  is  accomplished  by  assuming  that  the  large 

desorption  peak  at  approximately  150°K  is  composed  of  contributions  from 

nitrogen  bound  in  two  states,  with  different  values  of  E,  and  v, . We  know 

des  1 

the  shape  of  the  total  desorption  curve  peak  desorption 

temperature  (152°K)  from  Fig.  30.  All  that  remains  is  to  calculate 

the  values  of  the  desorption  parameters  of  the  low  temperature  state, 
designated  state  that  desorbs  at  a slightly  higher 

(2)  (2) 

temoerature  El  , v,  . The  values  of  the  desorption  parameters  are  obtained 
des  1 

by  varying  > ^des  ’ '^1^^’  until  the  best  fit  to  the 

observed  desorption  curve  is  obtained.  A copy  of  the  computer  program 
which  accomplishes  this  fitting  is  given  in  Appendix  B. 

The  result  of  the  separation  of  the  two  desorption  curves  is  shown  in 
Fig.  31.  Along  with  the  original  flash  desorption  curve  and  the  separated 
desorption  curves  is  the  sum  of  the  latter  two.  This  sum  does  not  add  up 
to  the  experimentally  observed  desorption  rate.  Since  there  are  a multitude 
of  different  crystal  planes  located  along  the  edge,  we  are  probably  dealing 
with  a number  of  different  states  that  vary  slightly  in  their  desorption 
energy  and  we  would  not  expect  this  simple  model  to  give  the  observed  curve. 
This  method  does,  however,  give  an  estimate  of  the  activation  energy  of 


desorption  from  the  (110),  which  we  find  to  be  Ej  = 5.4  + .76  kca 1/mole 
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Figure  31.  Decomposition  of  the  flash  desorption  curve  for  v nitrogen 

into  two  first  order  processes,  a)  =5.4  kcal/mole, 

(1)  _ .,10 


y"""  = 10^"  sec'"- , b)  = 8.7  kcal/mole,  = 1. 

^ 12  -1  des 

X 10  sec  , c)  the  sum  of  a and  b 
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with  = 10^*^  sec  for  a surface  coverage  of  0.32  monolayer.  The 

second  state  gives  = 8.7  + 1.9  kcal/mole  with  = 1.5  x 10^^  sec  ^ 

des  “ 1 

We  have  analyzed  the  desorption  rate  of  v nitrogen  from  the  (110)  plane  of 
tungsten  for  three  different  initial  concentrations,  0.45,  0.82  and  1.0 
monolayer.  The  results  of  this  analysis  are  shown  in  Table  1. 

Table  1 

Flash  Desorption  Parameters  of  v Nitrogen  from  the  W(llO)  Coverage 


Coverage  (monolayer)  (kcal/mole) 

E^^^  (kcal/mole 
des 

vj^^sec’ 

0.45  5.2 

8.8  X 10^ 

7.95 

2.7x  10^^ 

0.82  5.4 

1.0  X 10^° 

8.7 

12 

1.5  X 10 

1.0  5.1 

6.0  X 10^ 

3.05 

6.0  X 10^^ 

(2)  (2) 

The  values  for  ^nd  characteristic  of  the  higher  temperature 

state  are  very  similar  to  those  obtained  for  the  flash  desorption  of  v 
nitrogen  from  polycrystalline  tungsten  surfaces  [8],  It  may  be  that  the 
shoulder  we  observed  is  due  to  the  desorption  of  molecular  nitrogen  from  the 
edges  of  the  crystal.  The  major  (110)  peak  has  an  unusually  low  prefactor, 
when  compared  to  that  of  polycrystalline  and  other  planes  of  tungsten  [8,14]. 
The  activation  energy  and  prefactor,  however  are  both  very  similar  to  the 
only  work  that  has  been  done  with  v nitrogen  on  the  (110)  plane,  in  which 
the  activation  energy  for  desorption  was  found  to  be  E^^^  =6.4  kcal/mole 
for  = 10^^  sec  ^ [37]. 
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From  the  flash  desorption  curve  of  Fig.  30  we  also  notice  that  there 
is  some  3 nitrogen  formed;  the  desorption  temperature  for  this  state  is 
w 1500°K.  The  3 on  the  (110)  plane  amounts  to  about  7%  of  the  total  gas 
desorbed  as  determined  by  comparing  the  area  under  the  V and  3 parts  of  the 
desorption  curve.  Whether  the  formation  of  3 occurs  at  95°K,  or  whether 
it  occurs  while  the  crystal  is  being  flashed,  is  not  known.  From  the 
desorption  curves  little  information  can  be  obtained  on  the  desorption 
kinetics  of  3 nitrogen.  However,  we  do  notice  that  the  peak  desorption 

3 

temperature  of  3 nitrogen,  T^,  depends  upon  concentration,  decreasing  with 
increasing  coverage.  This  is  indicative  of  a second-order  reaction  [35,36]. 


IV. 2. B.  Structure  of  the  Gas  Laver 

When  a clean  tungsten  (110)  surface  (at  room  temperature)  is  exposed  to 
nitrogen,  also  at  room  temperature,  a distinctive  change  in  the  LEED  pattern 
is  observed.  Fig.  32  shows  the  LEED  pattern  as  the  crystal  is  exposed  to 
progressively  larger  amounts  of  nitrogen.  Prior  to  each  photograph  the 
crystal  is  heated  for  several  seconds  at  800°K  to  fully  develop  the  LEED 
pattern.  In  this  sequence  of  pictures  the  (0,0)  and  the  (1,1)  beam  are 
obscured  by  the  crystal  holder  and  one  of  the  heating  wires.  With  increasing 
coverage,  as  determined  by  Auger  spectroscopy,  and  heating  to  800°K  for 
several  seconds,  a p(2x2)  pattern  become.?  visible,  indicating  the  existence 
of  a new  surface  periodicity  with  a unit  cell  twice  that  of  the  old  surface 
net,  as  shown  in  Fig.  19.  The  basic  arrangement  is  one  in  which  the  nearest 
neighbor  sites  around  a given  nitrogen  atom  are  left  vacant.  It  should  be 
noted  that,  as  the  pattern  develops,  all  of  the  half  order  beams  split,  as 
shown  in  Fig.  32c,  d,  e,  and  f.  This  splitting  affords  some  insight  into 
the  growth  mechanism  of  a nitrogen  gas  layer. 


Changes  in  LEED  pattern  when  initially  clean  W(llO)  is  exposed  to 
progressively  greater  amounts  of  nitrogen.  In  Fig.  b,  c,  d,  and 
f,  the  sample  is  briefly  heated  to  800°K  for  several  seconds 
after  exposure.  The  pattern  in  e is  obtained  after  annealing  at 
room  temperature  for  15  minutes.  The  surface  concentrations  are, 
a)  clean,  b)  m 0.1  monolayer,  c)  «%5  0.18  monolayer,  d)  rt  0.2 
monolayer,  e)  rs  0.25  monolayer,  f)  0.25  monolayer. 


Figure  32 
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Beam  splitting  is  a well  known  phenomenon,  and  arises  fr-^m  the  forma- 
tion of  surface  domains  which  are  out  of  phase  with  each  other,  as  illustrated 
in  Fig.  19.  The  domains  of  Fig.  19  each  have  the  same  surface  structure, 
a p(2x2),  with  the  vector  d^  connecting  the  domains.  In  this  case  d^^ 
is  not  a net  vector  of  the  superlattice  and  the  domains  are  out  of  phase. 

From  the  discussion  in  Sec.  III. 2. A.,  the  only  beams  that  will  be  split  are 
those  which  satisfy  Eq . (18), 

hn  + km  = j (2i  +1)  i = integer,  (18) 

where  d^^  = na  + mb.  In  this  case,  the  only  LEED  beams  that  are  affected 
are  those  in  which  h and  k are  odd  multiples  of  j . When  the  antiphase 
boundaries  are  parallel  and  regularly  spaced,  the  beam  splitting  is  in  a 
direction  perpendicular  to  the  boundary.  A surface  structure  which  results 
in  the  observed  p(2x2)  LEED  pattern,  with  the  proper  beam  splitting  is 
illustrated  in  Fig.  19.  That  the  surface  structure  of  Fig.  19  does  indeed 
result  in  the  correct  LEED  pattern  was  verified  by  using  Eq . (15)  to  calcu- 
late the  beam  intensity  along  several  characteristic  directions.  These 
directions  are  shown  in  Fig.  32e.  The  surface  lattice  consisted  of  2 
p(2x2)  domains  separated  by  a translation  vector  d^  = ax,  as  illustrated 
in  Fig.  19,  each  domain  being  comprised  of  128  adatoms.  The  result  of 
this  calculation  is  shown  in  Figs.  33  and  34,  and  agree  with  the  LEED  pattern 
of  Fig.  32e.  The  location  of  the  LEED  beams  in  the  p(2x2)  structure  are 
listed  in  Table  2 . 

It  is  not  necessary  to  heat  the  crystal,  after  nitrogen  exposure,  in 
order  to  obtain  the  p(2x2).  Immediately  after  exposure  no  new  structure 
is  noticed  in  the  LEED  pattern;  however  after  10-20  minutes  at  room 
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Figure  33.  Calculated  scattered  intensity  along  line  (a)  of  Fig.  32e. 
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Figure  34.  Calculated  scattered  intensity  along  line  (b)  of  Fig.  32e. 
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Table  2 

LEED  Beam  Location 


P(2  X 2) 


(i,  j)  i = - 1,0,1 
j = - 1.0,1 


(n 

V2  ’ 2 i 


n = -l,l^,aj  'i,- 
, , ' douDlet 
m = - 1,1 


temperature,  a p(2x2)  pattern  is  observed  and  is  shown  in  Fig.  32e.  This 
behavior  is  understandable  if  initially  the  atoms  are  not  arranged  in  a 
periodic  structure  on  the  surface  but  become  ordered  on  the  surface  by 
diffusion.  An  estimate  of  the  root  mean  square  displacement  of  the  surface 
atom  based  on  these  LEED  measurements  will  be  given  in  Sec.  IV. 5. 

A slight  modification  of  the  standard  4-grid  LEED  system  (discussed 
in  Sec.  III.2.C.)  permits  us  to  measure  the  change  in  work  function 
accompanying  the  changes  in  the  LEED  pattern.  Fig.  35  shows  the  changes 
that  occur  in  the  characteristic  I -V  curve  as  the  surface  is  saturated 
with  3 nitrogen  at  room  temperature.  After  the  3 nitrogen  is  deposited,  a 
slight  change  in  the  I-V  curve  is  actually  noticed.  No  further  changes 
occur  after  the  crystal  is  warmed  and  a p(2  x2)  LEED  pattern  becomes  visible. 
The  initial  change  is  probably  not  due  to  a change  in  the  work  function, 
but  instead  arises  from  a slight  change  in  the  reflection  coefficient  r as 
discussed  in  Sec.  III.2.C.  Notice  that  the  various  curves  in  Fig.  35  are 
not  displaced  with  respect  to  one  another,  as  they  should  be  if  the  work 
function  changed,  but  actually  have  a slightly  different  slope 
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(in  the  linear  region).  This  is  just  the  effect  expected  if  the  reflection 
coefficient  changed.  At  lower  coverages  there  is  no  noticeable  difference 
in  any  of  the  I-V  curves.  This  agrees  with  earlier  studies  [1^,15]. in  which 
no  change  in  the  work  function  was  found  after  exposing  the  tungsten  (110) 
plane  to  nitrogen  at  room  temperature. 

IV. 2. C.  Desorption  Kinetics  of  Nitrogen  Adsorbed 
at  Room  Temperature 

Before  we  begin  to  study -the  surface  diffusion  of  nitrogen  on  the 
(110)  plane  of  tungsten,  it  is  first  necessary  to  characterize  the  surface 
gas  layer.  Not  only  is  it  important  to  know  which  state  the  gas  is  bound 
in,  but  it  is  also  of  considerable  interest  to  compare  the  binding  energy 
of  the  gas  on  the  surface  to  the  activation  energy  for  diffusion.  The 
most  common  method  of  obtaining  the  desorption  parameters  of  a gas  adsorbed 
on  a surface  is  to  heat  the  sample  at  a fixed  rate  and  measure  the 
accompanying  pressure  change,  either  with  an  ion  gauge  or  mass  spectrometer. 
This  was  just  illustrated  by  the  desorption  measurements  of  v nitrogen. 
Another  technique  that  can  be  used  to  follow  the  desorption  process 
involves  the  use  of  Auger  spectroscopy.  Basically  the  experimental  procedure 
is  simple.  A gas  layer  is  deposited  on  the  surface  and  .\uger  spectroscopy 
is  used  to  determine  the  amount  of  gas  remaining  on  the  surface  after 
the  crystal  is  heated  to  various  temperatures.  This  has  the  advantage  of 
determining  the  amount  of  gas  on  the  area  of  interest  only.  This  is  the 
first  time  that  Auger  spectroscopy  has  been  used  to  study  the  desorption 
kinetics  of  an  adsorbed  gas  layer.  To  determine  the  characteristics  of  the 
gas  layer  on  which  the  diffusion  studies  will  be  made,  we  compare  the 
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desorption  of  3 nitrogen,  using  Auger  spectroscopy  to  those  obtained  by 
flash  desorption  [13].  It  is  kno^m  that  nitrogen  adsorbed  at  room 
temperature  on  the  (110)  plane  of  tungsten  is  bound  in  the  9 state.  This 
state  exhibits  second-order  desorption  kinetics  with  a rate  given  by 
Eq.  (34),  so  that  the  concentration  C(t)  after  a heating  interval  of  t 
seconds  is  given  by 


l+C^v^t  expl,- 


'des  1 
RT  ( 
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where  is  the  initial  concentration  of  nitrogen  adsorbed  on  the  surface. 

Auger  spectroscopy  does  not  measure  the  surface  concentration  of 
nitrogen,  but  actually  gives  a signal  proportional  to  M(t),  the  number  of 
atoms  in  the  beam  area.  For  a uniform  surface  concentration  (within  the 
beam  area)  the  Auger  signal  is  proportional  to  C(t).  The  resultant  Auger 
signal  (proportional  to  M(t))  is  given  by 
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where  Mq  is  the  initial  number  of  atoms  sampled  and  M(t)  the  number  of  gas 
atoms  within  the  beam  after  heating  for  t seconds.  This  becomes 
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The  actual  experimental  procedure  is  as  follows.  After  the  surface 


is  cleaned  of  impurities,  the  crystal  is  exposed  to  a saturation  dose  of 
nitrogen  gas.  Since  the  sticking  coefficient  of  nitrogen  on  the  (110) 
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plane  at  room  temperature  is  < .004  , we  are  forced  to  expose  the  crystal 

“6 

;i  to  10  Torr  of  nitrogen  for  1000  seconds.  This  enables  us  to  obtain  a 

; reasonable  surface  concentration  that  can  be  detected  by  Auger  spectroscopy. 

: After  the  background  gas  is  removed,  an  Auger  scan  in  the  energy  range  of 

[ the  nitrogen  K , L.,  _ transition  is  made.  The  peak  to  peak  height  of 

1 

; this  Auger  transition  is  then  equal  to  oMg>  c being  a proportionality 

constant.  A number  of  Auger  scans  are  made  at  different  locations  on  the 
crystal  to  obtain  an  average  surface  concentration.  The  sample  is  heated 
to  a predetermined  temperature  for  a fixed  time  and  then  cooled  immediately; 

‘ another  Auger  scan  is  made,  to  obtain  M(t) . The  results  of  these  scans 

are  shown  in  Fig.  36.  This  procedure  is  repeated  a number  of  times  at 


different  temperatures.  Using  Eq . (40c)  we  see  that  a plot  of 

n / ^0 

Inj  ~ V w / s “1  j i . 1/T  should  yield  a straight  line  with  the  slope  proper- 
tional  to  the  activation  energy  of  desorption  and  the  y-intercept  equal  to 


Such  a plot  is  shown  in  Fig.  37.  We  find  the  activation  energy  of 

12 

desorption  (E,  ) to  be  76.9  +2.9  kcal/mole  and  C-v„  to  be  2.28  x 10 

^ des  — 0 2 

molecules  (cm'2  sec"^)  (5.6x10^^  < CQy2  9.3x10^^  molecules  (cm  ^ sec  ^) . 

The  method  used  to  determine  the  errors  in  and  E,  is  discussed  in 

.Appendix  C.  These  values  are  in  good  agreement  with  the  values  obtained  by 

Tamm  and  Schmidt  [13]  who,  using  flash  desorption  techniques,  found  E,  = 79 

des 

kcal/mole.  When  the  desorption  results  are  plotted  as  la'"^  ln(MQ/M)  ' vs.  1/T 

(i.e.  n=l)  the  resultant  cur'/e  has  a greater  standard  deviation  associated 

with  it  than  for  the  n=2  case.  The  concentration  dependence  of  the  peak 

3 

desorption  temperature  T^  in  Fig.  30  is  also  indicative  of  a second  order 
process.  Although  our  results  are  far  from  conclusive,  it  appears  that 


Auger  Intensity  (Arbitrary  Units) 
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Figure  36.  a)  Variation  in  concentration  across  nitrogen  layer  after  an 

exposure  of  10”^  Torr  for  1000  seconds  at  room  temperature  FI . 
b)  Variation  in  nitrogen  concentration  after  heating  (a)  to 
1048°K  for  1000  seconds,  A. 
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nitrogen  on  the  (110)  plane  of  tungsten  displays  second-order  desorption 
kinetics  . 


IV.3.  Nitrogen  Lavers  Produced  bv  Electron  Bombardment 
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IV. 3. A.  Electron  Stipulated  Dissociation 

Whenever  a surface  is  probed,  be  it  with  electrons,  atoms  or  photons, 
there  is  always  the  possibility  that  the  method  of  examination  will 
disturb  what  is  being  sampled.  This  has  been  found  to  be  the  case  when 
electrons  interact  with  a number  of  surface  gas  layers  [37-41].  We  actually 
use  the  interaction  of  an  electron  beam  with  an  adsorbate  to  overcome  one 
of  our  basic  problems,  that  of  creating  a localized  deposit  of  nitrogen  on 
the  surface.  The  interaction  of  an  electron  beam  with  nitrogen  bound  in 
the  Y state  appears  to  cause  dissociation  of  the  molecular  nitrogen. 

A possible  mechanism  for  this  dissociation  is  represented  schematically  in 
Fig.  38.  Electron  impact  excites  the  nitrogen  molecule  into  an  anti- 
bonding state  in  which  it  dissociates  on  the  surface.  The  component  atoms 
then  move  away  from  one  another  and  are  bound  as  atoms  on  the  surface. 

The  technique  used  to  create  these  localized  deposits  is  illustrated  in 
Fig.  39.  Initially  the  sample  is  cooled  to  95°I<  and  exposed  to  molecular 
nitrogen.  At  this  temperature  the  nitrogen  is  adsorbed  in  the  ^ state. 

The  sample  is  then  irradiated  with  an  electron  beam,  from  the  glancing 
incidence  Auger  g’un,  at  an  energy  of  2.5  KeV  and  a current  of  100  -A  in  a 
beam  of  0.1  cm  diameter.  Dissociation  of  the  molecular  nitrogen  occurs 
wherever  the  electron  beam  strikes  the  surface.  The  crystal  is  then 
warmed  above  « 120°K  and  the  weakly  bound  v nitrogen  desorbs,  leaving  only 
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nitrogen  on  the  (110)  plane  of  tungsten  displays  second-order  desorption 
kinetics . 

IV .3.  Nitrogen  Lavers  Produced  bv  Electron  Bombardment 
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incidence  Auger  g’un,  at  an  energy  of  2.5  KeV  and  a current  of  100  in  a 
beam  of  0.1  cm  diameter.  Dissociation  of  the  molecular  nitrogen  occurs 
wherever  the  electron  beam  strikes  the  surface.  The  crystal  is  then 
warmed  above  « 120°K  and  the  weakly  bound  Y nitrogen  desorbs,  leaving  only 
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Figure  39.  Operations  to  form  localized  deposits,  a)  Crystal,  at  95°K  is 
exposed  to  nitrogen,  forming  a layer  primarily  of  v nitrogen; 
b)  electron  irradiation  along  strip  enhances  concentration  of 
nitrogen  atoms,  c)  warming  to  T > ’l20OK  removes  v nitrogen, 
leaving  the  tightly  bound  atomic  nitrogen. 
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the  more  tightly  bound  atomic  nitrogen.  A high  energy  primary  beam  is  not 
necessary  to  form  these  deposits.  Using  the  glancing  incidence  Auger  gun 
we  have  formed  deposits  with  beam  energies  lower  than  1 KeV,  and  by  using 
the  LEED  electron  gun  deposits  have  been  formed  at  100  volts  with  1 jA 
beam  current. 

It  should  be  kept  in  mind  that  this  method  of  creating  surface  deposits 
is  only  practical  if  two  conditions  are  met.  First,  the  molecular  gas  must 
not  dissociate  spontaneously  on  the  surface  (or  at  worst  it  must  dissociate 
slowly),  and  secondly  it  must  have  a small  cross  section  for  electron 
stimulated  desorption  and  a high  cross  section  for  electron  stimulated 
dissociation . 

Fig.  40  shows  a secondary  electron  emission  photograph  of  a line  source 
of  atomic  nitrogen  on  the  (110)  plane,  at  95°K.  This  deposit  is  created  by 
following  the  procedure  described  earlier  in  this  section.  Adjacent  to  the 
line  of  atomic  nitrogen  is  the  undissociated  molecularly  bound  y nitrogen. 

If  the  crystal  is  warmed  above  120°K,  the  y desorbs  and  all  contrast  between 
the  line  of  atomic  nitrogen  and  the  surrounding  is  lost.  It  should  be  kept 
in  mind  that  while  this  picture  is  being  taken  the  y nitrogen  is  undergoing 
conversion  to  the  atomically  bound  state,  due  to  continual  irradiation 
from  the  Auger  gun.  Since  the  emission  current  necessary  to  produce  this 
picture  is  only  10  uA  and  the  beam  is  rastered  over  a large  area,  3-4  cm", 
the  effective  current  density  is  sufficiently  small  that  we  were  unable 
to  detect  any  y to  3 conversion  occurring  while  viewing  the  sample.  We  are 
not  restricted  to  line  source  as  the  starting  geometry.  Fig.  41  shows  a 
point  source,  approximately  0.1  cm  in  diameter;  also  in  this  picture  are 
the  .005"  thermocouple  wires  and  the  crystal  holder.  Larger  deposits  can 


Secondary  electron  image  of  a line  of  atomic  nitrogen  deposited 
on  the  (110)  plane  of  tungsten  at  95°K.  Primary  electron  beam 
current  is  10  m-A  at  2500  volts.  The  crystal  was  not  heated  to 
remove  the  y nitrogen  and  the  line  is  surrounded  by  y nitrogen. 


Figure  41,  Secondary  electron  image  of  a point  source  (~  .1  cm  in  diameter) 
of  atomic  nitrogen  on  the  (110)  plane.  The  crystal  is  cut  in 
an  elliptical  shape  with  a minor  axis  diameter  of  0,25". 


86 


be  created,  as  shown  in  Fig.  42,  by  changing  the  raster  pattern  of  the 
Auger  gun  during  the  irradiation  stage.  In  the  actual  experiment  once  a 
deposit,  such  as  shown  in  Fig.  40  and  Fig.  42,  is  created,  the  temperature 
is  quickly  raised  to  » 300°K.  Thermal  equilibration  takes  approximately 
2-3  minutes,  and  leaves  a localized  surface  deposit  of  atomic  nitrogen. 

The  amount  of  chemisorbed  nitrogen  deposited  by  electron  stimulated 
dissociation  can  be  controlled  by  one  of  two  methods.  The  first  is  by 
depositing  vatq/ing  amounts  of  v nitrogen  on  the  surface  (the  concentration 
being  less  than  one  monolayer)  which  is  then  bombarded  to  give  complete 
conversion  to  the  atomically  bound  state.  In  Sec.  IV. 2. A.  we  found  that 
when  y nitrogen  is  adsorbed  on  the  surface  there  is  some  conversion  into 
the  atomic  ® state.  At  a saturation  coverage  of  v nitrogen  this  conversion 

t 

amounts  to  ~ 0.17  monolayer.  Therefore,  by  using  the  method  described 
above,  there  will  always  be  some  3 nitrogen  on  the  surface.  The  concen- 
tration of  this  3 layer  will  be  dependent  upon  the  amount  of  v nitrogen 
originally  deposited.  This  method  has  the  advantage  that  the  ratio  of 
atomic  nitrogen  outside  the  irradiated  area  to  that  within  the  irradiated 
region  will  be  minimal,  » 30%. 

The  other  and  more  reproducible  method  is  to  deposit  a saturation 
amount  of  y nitrogen  and  then  control  the  current  density  to  the  v layer; 
this  then  controls  the  total  amount  of  nitrogen  dissociated  by  controlling 
the  number  of  electrons  which  strike  the  surface.  The  second  method  is 
the  one  we  used  in  forming  our  localized  deposits.  Its  major  drawback 
is  that  the  amount  of  atomic  nitr-igen  outside  the  irradiated  area  will 
always  be  determined  by  the  original  amount  of  y deposited,  irrespective 
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Figure  42.  Secondary  electron  image  of  a rectangle  of  atomic  nitrogen 
deposited  on  the  (110)  plane  of  tungsten. 
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of  the  amount  of  irradiation.  However,  this  method  is  capable  of  giving 
greater  reproducibility  in  the  deposit  concentration.  The  only  drawback 
being  a higher  concentration  ratio  of  atomic  nitrogen  outside  the  line 
to  the  inside  the  line. 

Auger  spectroscopy  is  also  used  to  compare  the  surface  coverage  for 
nitrogen  adsorbed  at  room  temperature  and  for  nitrogen  deposited  by 
eletron  stimulated  dissociation.  The  maximum  surface  concentration  of 
atomic  nitrogen  after  irradiation  is  » 2.2  times  the  room  temperature 
saturation  value.  To  determine  if  this  truly  is  the  maximum  attainable 
concentration,  we  deposit  a square  of  atomic  nitrogen,  redose  the 
crystal  and  again  irradiated  the  same  area.  An  Auger  spectmim  taken  of  the 
irradiated  area  shows  no  noticeable  increase  in  the  Auger  signal.  It 
therefore  appears  that  the  saturation  coverage  of  electron  dissociated 
nitrogen  is  indeed  2.2  times  the  saturation  coverage  at  room  temperature. 

Auger  spectroscopy  by  itself  is  unable  to  give  absolute  surface 
coverages.  However,  by  using  the  LEED  results  of  Sec.  IV. 2. B-  in  conjunc- 
tion with  Auger  spectroscopy,  an  estimate  of  the  surface  coverage  can  be 
made.  From  the  structure  shown  in  Fig.  19  we  estimate  the  surface  coverage 
after  the  crystal  is  exposed  to  nitrogen  at  room  temperature  to  be  0.25 
monolayer.  Since  we  are  dealing  with  a number  of  surface  islands,  the 
average  surface  coverage  would  then  be  expected  to  be  slightly  less  than 
0.25  monolayer.  Tamm  and  Schmidt  [13],  using  flash  desorption  techniques, 
estimated  a surface  coverage  of  0.21  monolayer,  in  good  agreement  with  out 
LEED  results.  We  can  now  use  the  nitrogen  Auger  peak  corresponding  to 
saturation  of  the  room  temperature  crystal  as  a reference  mark,  equal  to 
0.25  monolayer.  The  saturation  coverage  of  the  electron  irradiated  nitrogen 
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14  2 

is  thus  0.55  monolayer  (7.7x10  a toms /cm  ). 

IV. 3. B.  Cross  Section  for  Electron  Stimulated  Dissociation 
A knowledge  of  the  efficiency  with  which  Y nitrogen  is  dissociated 
by  electron  irradiation  should  be  helpful  in  understanding  the  mechanism 
of  electron  stimulated  dissociation.  However,  estimating  the  dissociation 
cross  section  is  difficult.  This  is  due  to  the  many  secondary  electrons 
emitted  from  the  surface  which  can,  and  probably  do,  dissociate  some  of 
the  surface  y nitrogen.  To  obtain  a crude  estimate  of  the  dissociation 
cross  section  we  neglect  the  effect  of  the  secondaries  entirely.  Under 
these  conditions  the  rate  of  dissociation  of  y nitrogen  is  given  by 


dM*  ^ ^e^'^^diss 
dt  e 


(41) 


where  je  is  the  current  density,  e the  charge  of  the  electron,  M'  the 

number  of  molecules  on  the  surface,  and  Q,.  the  dissociation  cross 

^diss 

section.  Solving  this  equation,  we  find  that  the  number  of  molecules 

I 

remaining  undissociated  on  the  surface  after  time  t ia 


M' 
u _ 

vT 

Mq  L e 


' ^e^diss^'i 


(42) 


where  Mq  is  the  total  number  of  y molecules  on  the  surface  at  time  t =■  0 
that  can  undergo  dissociation. 


On  irradiating  a . 1 cm  area  saturated  with  Y nitrogen  at  100  uA  for 
”u 

100  seconds  we  find  — = 0.25.  Use  of  Eq.  (42)  gives  a dissociation  cross 

^0 

• 18  2 

section  of  Q,.  **3x10  cm^ . This  compares  quite  well  with  the  work  of 

diss 
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Madey  and  Yates  [38]  and  Emrich  [41]  who  found  ^ 10  cm  . 

IV. 3. C.  Electron  Stimulated  Desorption 

During  both  the  diffusion  and  desorption  measurements  we  are  constantly 
probing  the  surface  with  an  electron  beam.  It  is  therefore  important  to 
determine  if  the  interaction  of  the  electron  beam  with  the  atomically  bound 
nitrogen  results  in  the  desorption  of  the  gas  layer.  To  determine  the 
effect  of  an  electron  beam  on  the  adlayer  we  create  a large  deposit  of  B 
nitrogen,  such  as  shown  in  Fig.  42.  This  deposit  is  then  irradiated,  at  a 
point  in  the  center  after  warming  to  w 300Ok,  with  a 2500  eV  beam  of  electrons 
at  150  uA  for  100  minutes.  Fig.  43  shows  the  variation  in  Auger  intensity 
after  extended  electron  bombardment.  After  100  minutes  of  irradiation  we 
find  no  change  in  the  concentration  of  the  adlayer. 

If  adsorption  from  the  gas  phase  during  irradiation  is  negligible, 
then  the  desorption  rate  is  given  by 


dM' 

dt 


VW*' 


(43) 


where  M'  is  the  surface  concentration  at  time  t and  Q^^g  the  desorption 

cross  section.  Solving  Eq.  (43)  for  Q^gg  (assuming  that  Auger  spectroscopy 

can  detect  a change  in  surface  concentration  of  .01  monolayer)  we  find 
-23  2 

Qdeg  <1,4x10  cm  . This  result  is  in  reasonable  agreement  with  the 

results  of  Emrich  [35,41]  and  Madey  and  Yates  [38]  who  found 
*2^  2 

< 5 X 10  cm^  for  B nitrogen  adsorbed  on  polycrystalline  tungsten. 


Auger  Intensity 
(Arbitrary  Units) 
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IV  .3.0.  Work  Function  and  LEED  of  Electron  Dissociated  Nitrogen 

By  using  electron  stimulated  dissociation  we  are  able  to  follow  the 
changes  that  occur  in  the  LEED  pattern  of  a clean  tungsten  (110)  surface 
as  the  surface  concentration  of  nitrogen  increases  from  0 to  0.55  monolayer. 
To  achieve  this,  the  crystal  is  cooled  to  95°K  and  electron  irradiated  to 
form  a localized  deposit.  After  a thermal  anneal  to  800°K  for  several 
seconds,  a change  in  the  clean  tungsten  (110)  LEED  pattern  is  observed. 

Heating  to  800°K  does  not  result  in  the  desorption  of  3 nitrogen. 

This  is  determined  by  comparing  the  Auger  peak-to-peak  height  of  the 
nitrogen  K L2  ^ ^2  3 ^^^^sition,  before  and  after  the  sample  is  heated.  In 
both  cases  the  Auger  intensity  is  found  to  be  identical.  For  surface 
concentrations  less  than  0.25  monolayer,  which  is  the  saturation  coverage 
of  the  (110)  plane  at  room  temperature,  the  p(2x2)  pattern  of  Section 
IV. 2. B.  is  found.  For  surface  coverages  higher  than  this,  another  structure, 
called  structure  A,  appears.  Figure  44  shows  the  changes  that  occur  in 
the  LEED  pattern  for  progressively  higher  surface  coverages.  Heating  the 
crystal  is  not  necessary  to  develop  structure  A.  After  15  to  20  minutes 
at  room  temperature,  structure  A became  visible.  Very  rarely,  1 out  of 
15  times,  a (2x4)  pattern  appeared.  Although  we  do  not  have  sufficient 
data  to  uniquely  determine  the  atomic  arrangement  within  the  unit  cell  of 
either  structure,  possible  repeating  units  are  shown  in  Fig.  45  and  Fig.  46 
for  structure  A and  the  (2x4)  pattern,  respectively.  The  domains  in 
Fig.  46  are  related  by  a rotation  and  not  a translation.  Such  domains  do 
not  cause  splitting  of  any  LEED  beams.  Only  domains  of  the  same  rotational 
type  will  interfere  with  each  other.  We  have  verified  this  by  the  computer 
program  given  in  Appendix  B,  using  the  surface  structure  shown  in  Fig.  46. 


Figure  44,  Changes  in  the  LEED  pattern  of  an  initially  clean  tungsten  (110) 
plane  after  deposition  of  an  adlayer  of  6 nitrogen,  followed  by 
increasing  periods  of  electron  irradiation.  With  the  exception 
of  Fig,  36d,  the  crystal  was  heated  to  800°K  for  several  seconds 
to  fully  develop  the  LEED  pattern.  Structure  A is  shown  in  d 
and  e,  along  with  the  p(2  x 4)  in  b and  c.  The  surface  concen- 
trations in  monolayers  are  a)  clean,  b)  ss  0,14,  c)  - 0.25, 
d)  0,55,  e)  ^.5  0,55,  f)  a (2  x 4)  pattern  seen  occasionally 
in  place  of  d and  e,  with  a surface  coverage  of  0.55, 
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Figure  45.  A possible  unit  cell  capable  of  producing  structure  A. 
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The  resultant  intensity  pattern  is  shown  in  Fig.  47,  along  the  direction 
shown  in  Fig.  44f.  We  find  that  no  splitting  results  and  that  the 
resultant  LEED  pattern  is  the  sum  of  the  LEED  patterns  from  each  domain. 

We  also  calculate  the  LEED  pattern  for  the  surface  structure  illustrated 
in  Fig.  45.  The  patterns  along  three  directions  in  reciprocal  space, 
labeled  a,  b and  c in  Fig  . 44e,  are  calculated  and  shown  in  Figs.  48,  49 
and  50,  respectively.  These  calculations  were  performed  for  a two- 
dimensional  grid  consisting  of  over  700  adatoms.  All  of  the  LEED  beams 
of  Fig.  44e  and  Fig.  44f  are  accounted  for.  However,  it  should  be  noted 
that  the  computed  intensities  vary  over  two  orders  of  magnitude,  which 
surely  is  not  the  case  in  the  experiments  shown  in  Fig.  44e. 

The  relative  intensities  of  the  calculated  LEED  beams  can  be  changed, 
while  maintaining  the  same  LEED  pattern,  by  changing  the  arrangement  within 
the  unit  cell  of  Fig.  45.  Let  us  use  the  4 atoms  in  the  comer  of  the 
parallelogram  in  Fig.  45  as  our  new  unit  cell.  A calculated  LEED  pattern 
similar  to  Figs.  48-50  is  obtained,  except  that  all  LEED  intensities  are 
now  identical.  Our  calculations  only  provide  information  on  the  positions 
of  the  LEED  beams.  To  account  for  the  intensities,  a proper  multiple 
scattering  calculation  is  necessary.  Table  3 provides  a listing  of  the 
location  of  all  LEED  beams  in  both  Structure  A and  the  (2  x4)  pattern. 

In  Sec.  IV2.B.  we  found  the  adsorption  of  P nitrogen  at  room  tempera- 
ture did  not  change  the  work  function  of  the  clean  tungsten  surface.  To 
determine  the  effect  of  electron  dissociated  nitrogen  on  the  clean  tungsten 
work  function  we  again  use  the  retarding  diode  analyzer.  Fig.  51  shows  the 
characteristic  I-V  curves  for  a clean  tungsten  surface,  for  a surface 
saturated  with  nitrogen  by  electron  stimulated  dissociation,  and  after 
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Figure  50.  Calculated  scattered  intensity  along  line  (c)  of  Fig.  44e. 
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Figure  51.  I-V  curves  for  a)  an  initially  clean  tungsten  (110)  surface 
b)  the  (110)  plane  with  a saturation  coverage  of  atomic 
nitrogen  0.55  monolayer  □ ; c)  (b)  flashed  to  800°K  for 
several  seconds,  structure  A visable  after  the  flash,  0. 


LEED  Beam  Location 


102 


103 


heating  to  800°K  for  several  seconds,  also  at  saturation.  There  is  no 
difference  in  the  clean  and  nitrogen  covered  I-V  curves  and  a slight  shift 
after  the  crystal  is  heated.  After  heating,  the  clean  tungsten  LEED 
pattern  changes  to  that  of  structure  A.  The  change  in  work  function  that 
occurs  is  less  than  - .027  eV,  and  in  fact  the  shift  revealed  in  these 
experiments  is  probably  due  to  a change  in  the  reflection  coefficient  after 
the  nitrogen  becomes  ordered  on  the  surface,  rather  than  to  a change  in 
the  work  function.  This  appears  to  be  the  case  since  the  various  I-V 
curves  are  not  parallel  to  one  another,  but  instead  seem  to  have  a different 
slope,  which  is  characteristic  of  a change  in  the  reflection  coefficient. 

IV. 3. E.  Thermal  Desorption  of  Electron  Dissociated  Nitrogen 

Throughout  the  earlier  sections  we  have  tacitly  assumed  that  nitrogen 
formed  on  the  surface  by  electron  bombardment  is  bound  in  the  same  state 
as  nitrogen  adsorbed  at  room  temperature.  However,  on  other  crystal 
planes,  for  instance  the  (100)  and  on  polycrystalline  samples,  this  is  not 
the  case  [37,39].  To  verify  that  the  nitrogen  on  the  surface  after  electron 
irradiation  is  actually  in  the  3 state,  we  compare  the  desorption  kinetics 
of  3 nitrogen,  formed  by  adsorption  at  300°K,  as  described  in  Sec.  IV. 2. C., 
to  that  of  the  electron  irradiated  nitrogen.  To  do  this  we  create  a large 
deposit  of  atomic  nitrogen  by  irradiation,  as  shown  in  Fig.  42,  and  follow 
the  desorption  from  the  center  of  the  deposit  as  a function  of  heating 
time  and  temperature,  using  Auger  spectroscopy.  The  deposit  is  sufficiently 
large  that  the  central  portion  is  unaffected  by  diffusion  from  the  edges 
of  the  deposit. 
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It  is  first  necessary  to  determine  the  order  of  desorption,  n.  This 

can  be  done  knowing  the  initial  concentration  and  the  time  tj^  needed  to 

2 

desorb  one-half  the  initial  surface  concentration.  The  desorption  process 
can  be  expressed  in  the  form 


= k(MQ -M)''  , (44) 

where  k is  a constant,  M the  surface  coverage  and  n the  order.  Upon 
integration  of  Eq.  (44)  and  substitution  of  M = obtain  for  a 

first  order  process 


and  for  n 1 


^ (2^'^-l) 

^ k(n-l)MQ"‘^ 


(45) 


(46) 


In  general,  the  time  required  for  one  half  the  original  concentration  to 
desorb  is  given  by 


ti 


f (n,k) 


(47) 


By  applying  Eq.  (46)  to  two  pairs  of  t^  and  M(t)  data  points,  an  expression 
for  the  desorption  order  can  be  obtained,  namely 


n 


log  t ^ - log  1 1 
log  Mq  - log  Mq 


(48) 
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where  tl  is  the  time  required  to  reduce  the  surface  concentration,  to 
2'  0 

one  half  its  initial  value  at  a fixed  temperature.  By  depositing  two 

different  concentrations  (Mq  and  Mq)  and  determining  the  time  needed  to 

remove  one-half  the  initial  concentration  the  order  can  be  calculated. 

Figure  52  shows  the  results  of  3 different  experiments  at  a temperature  of 

1072°K,  used  to  determine  t^.  Table  4 gives  the  values  of  the  initial 

2 

surface  concentration  (normalized  as  described  earlier)  and  the  time  needed 
to  reduce  the  surface  concentration  to  j its  initial  value.  Using  Eq . (48) 
the  desorption  order  is  2.2  + .1. 

The  measurement  suggests  that  desorption  is  indeed  second  order.  The 
activation  energy  for  desorption  can  now  be  found,  by  following  the  desorp- 
tion process  as  a function  of  time  and  temperature.  Figure  53  shows  the 
concentration  variation  across  the  central  portion  of  the  square  deposit, 
initially  at  saturation  and  then  after  heating  for  300  seconds  at  1048°K 
and  100  seconds  at  1095°K.  This  is  repeated  a number  of  times  at  different 

r 1 / ^0  ''1 

temperatures  and  for  different  times.  When  plotted  as  In  [ J v "m”  ''^s. 

1/T,  a straight  line  results.  This  is  shown  in  the  upper  portion  of  Fig.  37. 

The  slope  of  this  line  is  proportional  to  the  activation  energy  of 

desorption,  the  y intercept  to  CqU2’  activation  energy  to  be 

13 

79.2  + 5.1  kcal/mole  and  the  pre-exponental  constant  to  be  4.63  x 10 
-2  -1 

molecules  (cm  sec  ).  Also  plotted  in  Fig.  37  are  the  desorption  results 
for  nitrogen  adsorbed  at  room  temperature,  regular  3 nitrogen.  The  two 
curves  have  similar  slopes,  but  differeny  y-intercepts . The  difference 
in  the  intercepts  is  due  in  part  to  the  differences  in  initial  concentrations. 
The  saturation  coverage  of  electron  dissociated  nitrogen  is  2.2  times  that 
of  nitrogen  adsorbed  at  room  temperature.  A possible  reason  for  the 


Heating  Time  (sec) 


AP-828 


Figure  52.  Plot  used  to  determine  t^.  The  initial  concentrations  were 

a)  0.55  monolayer  A,  b)  « 0.52  monolayer  □ , c)  « 0.41 

monolayer  0,  with  ti  equal  to  316,  302,  and  458  seconds 
respectively.  ^ 
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Table  4 

Time  Necessary  to  Desorb  Half  the 

Initial  Concentration 

M/Mji 

t^  (sec) 

1 + .07 

316  + 16 

0.95  + .07 

302  + 22 

0.74  + .14 

458  + 14 

Auger  Intensity  (Arbitrary  Units) 
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Figure  53.  a)  Variation  in  surface  concentration  of  electron  dissociated 
nitrogen  across  the  crystal  face  Q . b)  Variation  in  surface 
concentration  after  heating  to  1048OK  for  300  seconds  4. 
c)  Variation  in  surface  concentration  after  haating  to  1095OK 
for  100  seconds  0. 
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difference  in  saturation  coverage  between  electron  dissociated  nitrogen 
and  regular  3 nitrogen  is  that  the  adsorption  process  at  room  temperature 
becomes  exceedingly  slow  for  coverages  > 0.25  monolayer.  Electron 
stimulated  dissociation  replaces  the  normal  adsorption  mechanism  with  a 
more  efficient  one  and  we  obtain  a higher  surface  concentration. 

We  have  found  that  the  nitrogen  which  is  deposited  on  the  surface  by 
electron  stimulated  dissociation  follows  the  same  desorption  kinetics 

13 

as  regular  3 nitrogen,  with  = 79.2  + 5.1  kcal/mole  and  4.63x10 

(4.14  X 10^^  < C^v^  < 5.19  X 10^^)  molecules  (cm  ^ sec  ^) . The  values  of 

12  "2  “1 

CqU^  regular  3 nitrogen  (2.28  x 10  molecules  (cm  sec  ))  and 
electron  dissociated  nitrogen  differ  by  a factor  of  9,  when  the  difference 
in  Cq  is  taken  into  account.  However,  the  differences  in  CqU^  well 
within  the  experimental  errors  associated  with  each  measurement.  From  this 
and  the  fact  that  the  I-V  curves  are  essentially  the  same  we  conclude  that 
the  electron  irradiated  nitrogen  is  also  bound  in  the  3 state. 

IV. 4.  Surface  Diffusion  Studies 

IV. 4. A.  Measurement  of  Diffusion  Profiles 

To  study  surface  diffusion,  a line  source  is  used  as  the  starting 
geometry,  as  it  is  both  easy  to  create  and  to  analyze.  After  the  v nitrogen 
is  deposited,  and  the  background  pressure  reduced  to  10  Torr,  an 
electron  beam  is  swept  across  the  crystal  face  for  300  seconds,  at  an 
emission  current  of  100  uA.  The  nondissociated  v nitrogen  is  removed  by 
warming  the  crystal  to  room  temperature.  A secondary  electron  image  of  a 
line  source,  produced  after  300  seconds  of  electron  bombardment  of  the  V 
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nitrogen,  is  shown  in  Fig.  40.  For  a concentration  independent  diffusion 
coefficient  and  a concentration  gradient  in  the  x direction,  we  know  that 


dC(x,t) 

dt 


D d C(x,t) 
dx2 


(49) 


where  D is  the  diffusion  coefficient  and  C(x,t)  the  concentration  after 
time  t at  a distance  x from  the  origin. 

For  a line  source  of  width  2h  the  concentration  C(x,t)  at  a distance  x 
perpendicular  to  the  center  of  the  initial  deposit  after  time  t is  [43] 


C(x,t) 


(50) 


where  Cg  is  the  initial  difference  in  concentration  inside  and  outside  the 
strip.  However,  since  we  are  sampling  a line  of  width  2h  with  a beam  of 
diameter  2h,  the  Auger  signal  must  be  intergrated  over  the  beam  dimension. 
Assuming  the  beam  is  uniform  and  has  a rectangular  cross  section,  with 
dimensions  b x 2h,  the  Auger  signal,  which  is  proportional  to  the  number 
M(x,t)  of  atoms  being  sampled,  is  given  by 


M(x,t)  _ jc_ 


M. 


4h 


erf 


( ^ ^ 

l^2h-x') 

, / 2h-x'\ 

^ \ 4h  y 

V ! 

+ erf 


2h_:^''i^ 


exp 


(51) 


where  Mg  is  the  number  of  nitrogen  atoms  initially  within  an  area  b x 2h,  and 
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M(x,t)  is  the  member  of  nitrogen  atoms  in  an  area,  centered  at  x,  of  size 
b X 2h  at  time  t. 

Thus  for  an  initial  deposit  in  which  the  concentration  is  a step 
function  of  the  distance  x,  the  Auger  signal  should  actually  appear 
triangular,  as  is  illustrated  in  Fig.  54.  When  diffusion  occurs,  the  base 
of  the  triangle  increases  as  material  diffuses  away  from  the  initial 
deposit  at  the  origin,  as  shown  in  Fig.  55. 

Once  a localized  deposit  (a  line)  of  nitrogen  is  formed  by  electron 
stimulated  dissociation,  scanning  Auger  spectroscopy  is  used  to  follow 
the  diffusion  process.  This  is  shown  schematically  in  the  lower  portion 
of  Fig.  54.  To  follow  the  diffusion  process,  a nxmber  of  Auger  spectra 
are  taken  at  discrete  locations  on  the  crystal  surface,  perpendicular  to 
the  line  of  nitrogen.  The  glancing  incidence  Auger  gun,  of  Fig.  54,  can 
position  the  electron  beam  anywhere  on  the  crystal  surface  by  the  use  of 
two  sets  of  electrostatic  deflection  plates.  Typically  a 2 mm  displacement 
in  the  beam  position  is  caused  by  5 volt  change  in  plate  potential.  At 
each  discrete  surface  position  an  Auger  scan  is  taken  within  the  energy 
range  370-390  V,  corresponding  to  the  K L-  -L,  - nitrogen  transition. 

In  this  manner,  an  integrated  concentration  profile  as  a function  of 
distance  from  the  center  of  the  line  is  obtained.  The  crystal  is  then 
resistively  heated  to  a predetermined  temperature  in  the  range 
800  < T < 900°K,  for  a fixed  period  of  time  (1000-3500  seconds  depending 
upon  the  temperature).  The  time  needed  to  reach  the  diffusion  temperature 
is  about  5 seconds,  and  amounts  to  .5%  of  the  shortest  diffusion  interval. 
After  the  predetermined  diffusion  time  has  elapsed,  the  crystal  heater 
is  turned  off  and  the  crystal  quickly  cools  to  room  temperature. 
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Figure  54.  a)  Schematic  representation  of  a scanning  Auger  microscope 

b)  Concentration  profile  perpendicular  to  a line  source. 

c)  The  observed  Auger  distribution. 
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Figure  55.  Expected  variation  in  Auger  intensity  perpendicular  to  initial 
line  source  of  atomic  nitrogen. 
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The  temperature  drops  to  400°K  in  about  20  seconds  and  within  15  minutes 
the  crystal  is  at  room  temperature.  At  that  moment  another  Auger  scan  is 
made  perpendicular  to  the  line  to  obtain  a new  concentration  profile. 

The  concentration  distribution  so  obtained  can  then  be  analyzed  to  obtain 
the  diffusion  parameters . 

After  measuring  a diffusion  profile,  the  crystal  is  cleaned  by 
flashing  to  a high  temperature,  a new  deposit  of  chemisorbed  nitrogen  is 
formed  by  electron  bombardment,  and  another  diffusion  sequence  is  initiated, 
but  this  time  at  a different  annealing  temperature.  During  all  of  the 
diffusion  measurements  the  time  and  temperature  of  the  diffusion  sequence 
were  chosen  to  insure  that  the  loss  of  material  through  desorption  was 
negligible . 

IV. 4. B.  Diffusion 

To  study  the  migration  of  nitrogen  atoms  over  the  tungsten  surface 
we  deposit  a line  of  atomic  nitrogen,  as  described  in  Section  IV. 3. A. 

Using  Auger  spectroscopy  we  are  able  to  follow  the  broadening  of  the  line 
as  diffusion  occurs.  Figure  55  shows  the  concentration  profiles  expected 
after  various  diffusion  intervals.  The  actual  spatial  variation  observed 
at  two  different  temperatures  is  shown  in  Figs.  56  and  57.  Plotted  in 
Fi'^.  56  is  the  concentration  perpendicular  to  the  initial  line  deposit 
before  and  after  the  crystal  is  heated  to  885°K  for  2000  seconds. 

Figure  57  shows  the  concentration  variation  across  a line  before  and  after 
heating  to  800°K  for  3500  seconds.  By  fitting  the  theoretical  profile. 


Eq . (51),  to  the  experimental  profile  we  are  able  to  obtain  a value 
of  the  diffusion  coefficient  D.  This  procedure  is  repeated  for  a number 
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Figure  56.  Spatial  variation  in  the  Auger  intensity  perpendicular  to  a 
line  of  atomic  nitrogen  of  width  .1  cm  for  t » 0 seconds  and 
after  heating  for  2000  seconds  at  885°K  A. 
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of  temperatures,  spanning  a range  of  100°K. 

Taking  the  diffusion  coefficient  to  have  the  usual  form 

- e’^-, 

D = Dq  exp'^  - , (52) 

with  the  activation  energy  of  diffusion,  and  Dq  a constant,  a plot 
of  ln(D)  vs.  1/T  should  result  in  a straight  line.  Such  an  Arrhenius  plot 
is  given  in  Fig.  58.  Also  shown  are  the  desorption  results  for  the 
electron  deposited  3 nitrogen.  This  is  the  first  time  desorption  and 


diffusion  measurements  have  been  made  on  the  same  crystal.  At  a satura- 
tion coverage  («  ,5  monolayer)  of  3 nitrogen  on  the  (110)  plane  of  tungsten; 

we  find  the  activation  energy  for  diffusion  E^  = 21  + 2.2  kcal/mole,  and 

m 

Dq  =4.66  X 10  cm^/sec  (1.25  x 10  ^ < Dq  < 1.7  x 10  ^ cm^/sec),  and  the 

13 

desorption  parameters  to  be  = 79.2  kcal/mole  and  CqV^  = 4.66  x 10 

(4.14  X 10^^  < CqV2  < 5.19  x 10^^)  molecules  (cm  ^ sec  ^),  as  previously 
reported.  The  method  by  which  we  determine  the  errors  in  both  the  diffusion 
and  desorption  parameters  is  discussed  in  Appendix  C. 

To  check  the  effect  of  concentration  on  the  diffusion  coefficient  we 
repeated  the  experiments  with  an  initial  coverage  only  3/4  of  saturation. 
This  reduced  surface  concentration  is  produced  by  exposing  the  cooled 


crystal  to  a saturation  dose  of  y nitrogen,  followed  by  sweeping  the 
electron  beam  across  the  crystal  face  for  100  seconds,  at  100  p.A  emission 
current,  instead  of  the  300  second  irradiation  in  the  previous  experiment. 


Since  we  started  with  a saturation  concentration  of  y nitrogen,  the  amount 

of  S nitrogen  adjacent  to  the  line  is  « 0.17  monolayer.  We  again  plotted 

ln(D)  vs.  1/T  and  the  results  are  shown  in  Fig.  59,  along  with  the 

-9  2 , 

saturation  coverage  results.  In  this  case  Dq  * 8.8  x 10  cm  /sec  and 
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of  surface  diffusion  0 and  of  desorption 
110).  Diffusion  parameters  found  are 
* 4.66  X 10"2  cm^/sec. 
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Figure  59.  Plot  of  ]n(D)  vs.  1/T  for  electron  dissociated  nitrogen  at 

a)  an  initial  surface  coverage  of  0.55  monolayer  A, 

b)  an  initial  surface  coverage  of  0.41  monolayer  0. 
Abscissa  at  left  for  high  concentration,  at  right  for  low 
concentration. 
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E = 21.27  +2.8  kcal/mole.  The  differences  in  the  two  measurements  are 
m ~ 

within  the  standard  deviations  of  each  of  the  activation  energies  and 
pre-exponentials.  The  diffusion  measurements  at  3/4  saturation  were  made 
with  a different  Auger  gun,  which  may  have  resulted  in  a slight  variation 
in  the  recorded  Auger  intensity.  To  check  the  effect  of  the  new  gun,  we 
repeated  some  of  the  diffusion  experiments  at  saturation  coverages,  and 
obtained  results  identical  to  those  with  the  old  gun.  Table  5 lists  the 
various  values  of  D obtained  for  different  heating  schedules  and  initial 
concentrations.  Although  our  range  of  concentrations  spanned  is  small,  it 
appears  that  diffusion  is  hardly  affected  by  changes  in  the  concentration 
of  the  nitrogen  adlayer. 

IV. 4. C.  Simultaneous  Diffusion  and  Desorption 

We  have  shown  that  electron  bombarded  nitrogen  desorbs  with  the  same 
kinetics  as  does  P nitrogen.  However,  flash  desorption  only  gives 
information  about  the  surface  species  at  the  desorption  temperatures,  about 
1100°K  in  this  case.  It  could  still  be  argued  that  at  the  lower  temperatures 
at  which  diffusion  occurs,  900°K,  the  electron  deposited  nitrogen  is  in 
still  another  binding  state.  It  could  be  that  only  at  1000°K,  just  prior 
to  evaporation,  does  it  convert  into  P and  desorb  as  such.  We  can  eliminate 
this  objection  by  examining  the  simultaneous  desorption  and  diffusion  of 
nitrogen  from  a line.  This  type  of  observation  is  bound  to  provide  infor- 
mation about  desorption  and  diffusion  of  nitrogen  in  the  same  state. 

Ideally,  it  is  desirable  to  operate  in  a temperature  and  time  regime  where 
the  mass  loss  from  the  initial  deposit  is  roughly  the  same  for  both 


mechanisms . 
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Table  5 

Diffusion  Parameters 


0.55  Monolayer  Coverage  0.41  Monolayer  Coverage 


0 

Temperature  ( K) 

2 

D(cm  /sec) 

Temperature  (°K) 

2 

D(cm  /sec) 

800 

7.56  X lO"® 

800 

1.169  X lO’^ 

815 

1.06  X lO'^ 

815 

1.6  X 10'^ 

835 

1.44  X lO'^ 

835 
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To  do  this,  an  appropriate  time  interval  and  temperature  must  first 
be  decided  for  the  experiment.  The  loss  of  material  by  diffusion  and 
second-order  desorption  is  given  by 


dC (x, t) 
dt 


D - kC^(x,t) 


dx 


(53) 


where  k = exp[E^gg/RT]  is  the  rate  constant  for  evaporation  and  C(x,t) 
the  surface  concentration.  This  equation  cannot  be  solved  exactly; 
however,  numerical  answers  can  be  obtained,  using  the  program  in  Appendix  B. 
The  calculated  concentration  is  then  integrated  over  the  beam  dimension  to 
obtain  which  is  proportional  to  the  Auger  signal.  A number  of 

different  heating  times  and  temperatures  were  tried  and  approximately  equal 
mass  loses  occured  for  t = 1000  seconds  and  T = 1000°K.  For  a heating  time 
t = 1000  seconds  and  temperature  T = 1000°K,  we  obtain  the  calculated  desorp- 
tion and  diffusion  curves  shown  in  Figs.  60  and  61.  These  figures  show  the 
variation  in  Auger  intensity  for  a 1 kcal/mole  change  in  the  activation 
energy  of  desorption.  In  this  temperature  and  time  range,  differences  in 
the  activation  energy  for  desorption  and  diffusion  should  be  readily 
discerned.  An  increase  of  one  kcal/mole  in  the  diffusion  energy  causes  a 
smaller  decrease  in  the  central  section  of  the  concentration  gradient  and 
a narrowing  of  the  base.  A similar  increase  in  the  desorption  energy 
causes  an  increase  in  the  overall  concentration  profile. 

The  result  of  an  actual  diffusion  and  desorption  anneal  at  T = 1000°K 
and  t = 1000  seconds  is  shown  in  Fig.  62.  The  best  fit  is  obtained  by 
using  the  diffusion  parameters  obtained  earlier  and  increasing  the  desorption 
energy  by  1.8  kcal/mole,  well  within  the  standard  deviation  of  the  desorption 
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Figure  60.  Calculated  changes  in  concentration  profile,  for  a 1 kcal/mole 
change  in  be  activation  energy  of  diffusion  when  diffusion 
and  desorption  occur  simultaneously.  E,  ■ 79.2  kcal/mole 
and  C..v_  = 4.63  x 10^^  /cm"^  sec”^^. 
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Figure  61.  Changes  in  the  concentration  profile,  when  diffusion  and 

desorption  occur  simultaneously,  resulting  from  a change  of 
1 kcal/mole  in  the  activation  energy  of  desorption;  = 21.0 
kcal/mole.  ^ 
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Figure  62.  Experimental  concentration  profile  for  simultaneous  diffusion 
and  desorption  from  a line  source. 
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energy.  The  diffusion  and  desorption  parameters  can  now  be  described  by 

^ -2  2 
E = 21.  kcal/mole,  = 4.66x10  cm“/sec,  E,  = 81  kcal/mole  and 
m 0 ’ des 

13  -2  -1 

= 4.63x10  atoms  (cm  sec  ).  We  now  have  diffusion  and  desorption 
occurring  simultaneously  and  the  activation  energies  are  almost  identical 
to  what  we  obtained  when  we  observed  each  event  occurring  separately. 

IV. 5.  Discussion 

Flash  desorption  and  Auger  spectroscopy  indicate  that  some  of  the 
molecularly  bound  nitrogen  undergoes  a conversion  to  the  more  tightly  bound 
3 state.  Whether  this  conversion  occurs  while  the  gas  is  on  the  surface 
at  the  lower  temperature  or  as  the  result  of  heating  to  some  higher 
temperature  is  unknown.  The  amount  of  3 found  in  this  manner  is  ^ 0.17 
monolayer.  There  is  the  possibility  that  lattice  steps  are  responsible 
for  this  V to  3 conversion.  From  our  LEED  measurements  we  know  that  the 
step  density  near  the  crystal  edge  is  « 1 step  for  every  30  surface  atoms, 
i.e.  the  terraces  are  30  atoms  wide.  Presumably  the  molecular  v nitrogen 
is  free  to  migrate  on  the  surface  at  95°K.  If  on  reaching  a step,  either 
by  impact  from  the  gas  or  by  diffusion  over  the  surface,  it  dissociates, 
then  we  would  have  two  atoms  bound  tightly  to  the  step,  and  no  further 
dissociation  could  occur  at  that  site.  The  maximum  concentration  in  this 
case  would  then  be  .07  of  a monolayer.  We  found  almost  twice  this  amount 
of  atomic  nitrogen  after  desorption  of  the  y nitrogen.  For  this  model  of 
dissociation  to  be  correct  every  step  atom  must  be  capable  of  dissociating 
two  molecules  of  nitrogen.  Another  possibility  is  that  the  step  density 


is  higher  than  1 in  30. 
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If  steps  are  necessary  for  the  dissociation  of  molecular  nitrogen  on 

the  (110)  plane,  then  at  room  temperature  the  gas  molecules  may  strike 

the  surface  and  remain  there  for  a short  period  of  time,  during  which  they 

can  diffuse  to  a step  and  dissociate.  If  during  their  residence  time 

they  do  not  migrate  to  a step,  they  leave  the  surface  as  a molecule.  This 

would  explain  the  low  sticking  coefficient  of  nitrogen  observed  on  the 

(110)  plane  at  room  temperature.  Once  a nitrogen  atom  is  on  a step,  it 

effectively  poisons  the  step  to  further  dissociation  and  only  by  diffusion 

from  the  step  can  another  molecule  dissociate  there.  The  relatively  low 

surface  coverage  observed  for  nitrogen  on  the  (110)  plane  at  room 

13  2 

temperature  (9.2x10  molecules/cm  ) [13]  would  then  be  the  result  of  a 
diffusion  limited  dissociation  process  at  the  steps.  By  electron  irradia- 
tion we  are  circumventing  the  normal  step  dissociation  mechansim  with 
another,  more  efficient  mechanism. 

As  mentioned  in  Sec.  IV. 2. 3 and  Sec.  IV. 3. D,  the  LZED  pattern  that 
develops  after  adsorption  of  nitrogen,  both  at  room  temperature  and  by 
electron  bombardment,  is  enhanced  by  a slight  thermal  anneal.  It  is  not 
necessary  to  heat  the  crystal  to  develop  either  the  p(2  x2)  or  structure  A. 
These  patterns  were  also  observed  after  the  crystal  was  allowed  to  remain 
at  room  temperature  for  10-20  minutes.  If  the  formation  of  the  I.ZZD 
pattern  is  a diffusion  controled  process,  an  estimate  of  the  root  mean 
square  displacement  necessary  to  form  the  LZZD  pattern  can  be  made. 
Initially,  the  nitrogen  atoms  are  assumed  to  be  arranged  randomly  on  the 
surface  and  during  the  10-20  minutes  needed  to  develop  a pattern  they 
diffuse  into  equilibrium  sites.  In  two  dimensions  the  root  mean  square 
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displacement  < r > of  an  atom  is  given  by 

< r > = v^^Dt  (54) 

By  extrapolating  the  high  temperature  diffusion  results  to  room  temperature 

-17  2 

we  obtain  a value  of  D = 2x10  cm  /sec.  This  gives  a root  mean  square 
displacement  of  a 10  atomic  spacings.  This  is  a somewhat  larger  displacement 
than  .might  have  been  expected,  considering  that  the  size  of  the  unit  cell 
of  the  adsorbate,  as  determined  by  LEED  measurements  is  only  a 4-8  atom 
spacings.  However,  the  experimental  error  associated  with  D is  sufficiently 
large  that  a root  mean  square  displacement  of  4 atoms  spacings  is  within 
the  experimental  error  of  the  measurements. 


SUMMARY 

We  have  developed  a method  of  creating  localized  deposits  of  atomic 
nitrogen  on  the  (110)  plane  of  tungsten  by  means  of  electron  stimulated 
dissociation.  Not  only  are  the  nitrogen  atoms  bound  in  the  3 state,  but 
with  a surface  concentration  2.2  times  that  which  can  normally  be  obtained 
by  room  temperature  saturation  of  the  (110)  plane.  The  nitrogen  which  is 
deposited  by  electron  irradiation  is  found  to  obey  second-order  desorption 
kinetics  and  to  desorb  with  an  activation  energy  of  79.2  kcal/mole.  This 
result  is  in  excellent  agreement  with  the  flash  desorption  measurements  for 
nitrogen  adsorbed  at  room  temperature. 

From  our  flash  desorption  studies  we  were  able  to  estimate  the  activa- 
tion energy  and  the  frequency  factor  for  desorption  of  > nitrogen  from  the 

(110)  plane  of  tungsten,  which  we  found  to  be  E,  ■ 5.4  + .76  kcal/mole 

aes  ~ 

and  V ■ 10  sec  Although  the  frequency  factor  is  unusually  low, 
when  compared  to  that  of  V nitrogen  on  the  other  planes  of  tungsten,  it 
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does  agree  very  well  with  the  only  other  work  which  has  been  done  on  the 
(110)  plane.  We  also  found  that  there  is  always  some  Y to  ^ nitrogen  at 

, room  temperature  is  w 0.25  monolayer. 

’ Our  LEED  observations  indicate  that  once  the  adsorbed  nitrogen  atoms 

equilibrate  on  the  surface  they  form  into  islands.  These  islands  have 
specific  orientations  and  the  borders  are  in  the  [lOO]  direction.  This  is 
not  unexpected,  since  island  growth  is  perhaps  the  most  common  form  of 
adsorption.  We  also  found  two  different  surface  structure.  One  corre- 
sponds to  a surface  concentration  of  < 0.25  monolayer  which  results  in  a 
p(2x2)  pattern.  This  pattern  forms  if  the  nearest  neighbor  sites  around 
any  one  nitrogen  atom  are  always  empty.  The  other  LEED  pattern,  structure 
A,  was  only  visible  for  surface  concentrations  > 0.25  monolayer.  We  also 
found,  on  occasion,  a structure  composed  of  rotated  (2x4)  domains.  ^ Both 
structure  A and  the  rotated  (2x4)  pattern  could  only  be  achieved  by 
electron  stimulated  dissociation. 

The  development  of  the  LEED  pattern  also  appears  to  be  a thermally 
activated  process.  After  the  nitrogen  is  deposited  on  the  surface  at  room 
temperature,  there  is  no  immediate  change  in  the  clean  tungsten  pattern. 

Only  after  the  crystal  is  heated  to  800°K  for  several  seconds  or  permitted 
to  remain  at  room  temperature  for  15-20  minutes  does  the  new  pattern  appear. 

m 

From  the  activation  energy  of  desorption  measured  in  our  experiments 
V we  can  deduce  the  binding  energy  of  a nitrogen  atom  on  the  surface.  The 

binding  energy  E^  is  given  by: 

■ r <X' 

> where  X'  is  the  dissociation  energy  of  the  molecule  in  free  space.  Using 

I 

I 


L 
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the  desorption  energy  found  earlier,  we  calculate  the  binding  energy  per 
atom  to  be  153  kcal/mole.  The  ratio  of  the  binding  energy  to  the  activation 
energy  for  surface  diffusion  is  7.4  to  1.  This  is  significantly  higher 
than  has  been  reported  on  rougher  surfaces.  From  field  emission  measure- 
ments, a ratio  of  approximately  4.5  to  1 appears  appropriate  for  nitrogen 
on  the  high  index  planes  of  tungsten  around  the  (110)  [4].  However, 
around  the  (111),  diffusion  occurred  with  the  same  activation  energy  as 
observed  here  on  the  (110).  These  comparisons  are,  at  best,  of  qualitative 
significance;  the  early  diffusion  measurements  were  only  semiqualitative , 
and  direct  observations  of  desorption  from  stepped  surfaces  are  not 
available.  It  is  of  interest,  however,  that  measurements  of  hydrogen 
diffusion  over  rough  surfaces  show  much  the  same  trend.  Values  of  » 7 
have  been  reported  for  the  binding  energy  to  the  diffusion  barrier  [3]. 

It  does  appear,  therefore,  that  in  diffusion  over  a smooth  surface  the 
activation  barrier  confronting  an  atom  is,  in  fact,  smaller  in  relation  to 
the  binding  energy  than  on  atomically  rough  surfaces. 

The  value  of  the  perfactor  Dq  in  the  surface  diffusion  coefficient 

for  nitrogen  atoms  on  the  tungsten  (110)  provides  some  information 

about  the  dynamics  of  atomic  motion  and  whether  surface  steps  are  limiting 

the  diffusion  process.  To  see  this  we  represent  the  surface  as  being 

L 

composed  of  a series  of  terraces  of  width  — , as  shown  in  Fig.  63.  For  a 
surface  atom  located  at  an  arbitrary  position,  i,  the  rate  constant 
describing  its  motion  to  the  right  is  The  atom  can  also  move  to  the 

left,  and  this  movement  is  assigned  the  rate  constant  u^,  as  indicated  in 
Fig.  63.  The  actual  values  of  the  rate  constants  depend  on  the  location 
of  the  surface  atom.  Jumps  on  the  flat  terraces  themselves  all  occur  at 
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1 Figure  63.  Schematic  drawing  of  a stepped  surface  with  adatom.  The 

‘ step  width  is  L/2  and  the  associated  rate  constants  of  the 

I adatom  are  a,  b,  c,  and  d. 
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the  same  rate,  independent  of  the  exact  position.  This  rate  is  denoted 
by  a.  At  the  edges,  hoveyer,  jump  rates  are  different,  and  are  assigned 
the  values  b,  c,  and  d shown  in  Fig.  63.  The  average  time  required  for 
an  atom  to  reach  either  of  two  end  points,  arbitrarily  chosen  as  x = 0 or 
X = L,  starting  at  x = i,  where  0 < i < L,  is  given  by  [44] 


where 


and 


Ai  + 


1 "^^2^2 


^2_+^l(l  ■*■22) 


^ m.ij=i'j  \ ’ 


(56) 


(57) 


L-1  m a 

B » E T7  ^ 

n=l  k=l 


(58) 


Using  the  rate  constants  of  Fig.  63,  for  an  atom  initially  at  point  i, 
where  L =■  2s  and  i > s,  we  obtain. 


Ai  - ^ (s-l)(2s-i-l)  +i  (2s-i-l) 


+ i U-l)  + ^ + i . 

c b 


(59) 


or 


12  2 
(2s  -6s44.i^+i 


(60) 


Specifically, 


(61) 


_ < 
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A similar  calculation  shows  that, 


(62) 


and 


Referring  to  Fig.  63,  we  see  that. 


"d  ^ "c  ■ "^a 


(63) 


Assuming  that  the  pre-exponential  factors  for  all  the  rate  expressions 
are  identical,  then. 


ad  = be 


(64) 


For  an  atom  initially  at  i = s,  it  then  follows  that 


B = (S-1+2-)  , 

s c 


(65) 


and 


[r  (-1) 


. s-1  . s^-3s+2 

- 

c 2a 


(66) 


The  transit  time,  now  follows  from  Eq.  (56)  as 


(67) 


To  visualize  the  effects  of  a step  which  acts  as  a trap  for  a migrating 
atom,  we  assume  that  a * c and  b ■ d;  the  transit  time  is  given  by 
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T « s (a+(s-l)b) 
s 2ab 


(68) 


where  a is  the  rate  constant  for  motion  on  the  terrace  and  b the  rate 
constant  for  removal  from  the  step  edge.  That  Eq.  (68)  is  reasonable  is 
clear  on  examining  a limiting  case,  which  corresponds  to  diffusion  over 
a flat  plane,  for  which  a * b.  In  this  case  Eq.  (68)  becomes 


s 2b  • 

and  the  mean  square  displacement  of  an  atom  in  time  t is  given  by. 


(69) 


(x^)  = -T-  (s2^)  = 2Dt  , 


(70) 


or 


D 


(71) 


which  is  just  the  result  obtained  for  a random  walk  experiment. 

For  the  case  of  interest  where  b « a,  diffusion  from  the  step  edge 
being  the  limiting  step,  Eq.  (68)  becomes. 


II 

T 

S 


(72) 


and 


— ^ (2s)^  = 2Dt  . (73) 

RT 

s 


The  diffusion  coefficient  is  now  given  by 


D 


bi^s 


U2^s  exp(-E’^/RT) 


f 


(74) 
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where  v is  the  frequency  factor  for  jumps.  The  pre-exponential  constant 
Dq  is  now  s.  That  is,  Dq  is  larger,  by  a factor  equal  to  the  terrace 
width  (s),  than  when  diffusion  occurs  over  a flat  surface.  In  our  experi- 
ment we  found  the  prefactor  to  be  quite  normal  and  of  a magnitude  expected 
for  an  atom  moving  over  sinusoidal  potential,  as  shown  in  Fig.  64,;  and 
we  therefore  conclude  that  trapping  at  lattice  steps  does  not  play  a 
dominant  role.  It  would  still  be  interesting  to  establish  the  possible 
contribution  of  steps  more  definitively,  both  by  measurements  of  diffusion 
along  and  perpendicular  to  the  step  edges,  and  by  artificially  producing 
steps  on  the  surface  during  the  initial  crystal  preparation. 

From  the  evidence  available  right  now,  in  any  event,  it  appears  that 
nitrogen  diffusion  on  the  (110)  plane  of  tungsten  occurs  in  quite  a simple 
fashion,  over  a rather  smaller  barrier  than  expected  from  studies  on  high 


index  surfaces. 
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APPENDIX  A 

Electron  Bombardment  Formation  of  Surface  Devices 

Electron  stimulated  dissociation  can  be  used  in  the  fabrication  of  a 
number  of  novel  surface  devices,  among  which  are  a surface  layer  mass 
memory  and  a possible  method  of  growing  localized  surface  deposits. 

Although  we  have  demonstrated  the  feasibility  of  depositing  nitrogen  on 
the  (110)  plane  of  tungsten,  it  may  be  possible  to  extend  this  technique 
to  other,  more  practical  systems. 

A surface  layer  mass  memory  can  be  fabricated  in  much  the  same  manner 

that  we  create  the  localized  nitrogen  deposits.  Sec.  IV. 3. A.  Now,  however, 

instead  of  depositing  a single  dot,  as  in  Fig.  41,  we  deposit  a binary 

code.  The  technique  by  which  this  code  is  deposited  is  shown  in  Fig.  A1 

where  a dot  represents  a "1"  and  a blank  space  represents  a "0". 

In  this  manner  a surface  memory  can  be  developed.  To  read  this  device, 

any  number  of  surface  sensitive  techniques  can  be  used,  such  as  Auger 

i n. .troscopy , secondary  electron  emission  or  electron  loss  spectroscopy, 

as  .liustrated  in  Fig.  A2.  An  actual  picture  of  such  a device  is  shown 

in  Fig.  A3,  which  was  taken  by  monitoring  the  secondary  electron  emission 

current.  A crude  estimate  of  the  time  necessary  to  write  on  this  device 

can  be  made  by  using  Eq.  (42)  and  the  value  of  the  dissociation  cross- 

• 18  2 

section,  Q..  , found  in  Sec.  IV. 3. B.,  which  is  Q,,  * 3 X 10  cm  . 

^diss  ^diss 

Assuming  a surface  concentration  of  0.28  monolayer  corresponds  to  a "1" 

and  that  the  electron  source  is  a field  emission  electron  gun,  capable  of 

•8  *9 

providing  a beam  current  of  10  -10  ampere  within  a beam  diameter  of 
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Figure  Al.  Writing  procedure  for  the  surface  layer  mass  memory,  a)  the 
molecular  gas  is  deposited  on  the  surface,  b)  followed  by 
local  irradiation  to  deposit  the  binary  code,  c)  warming  the 
crystal  to  remove  the  non-dissociated  molecules. 


Figure  A2 . Possible  techniques  (Auger  electron  spectroscopy,  electron 
loss  spectroscopy  and  seconday  electron  emission),  which 
may  be  used  to  read  the  surface  layer  mass  memory. 
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Figure  A3. 


Secondary  electron  emission  photograph  of  the  surface  layer  mass 
memory  device  with  a binary  code  on  the  surface. 


100  S,  the  writing  rate  is  ss  10  bits/sec.  Although  this  is  slow  by 

current  standards  this  technique  does  offer  the  possibility  of  achieving 

12  . 2 

very  high  bit  densities  10  bits/cm  , considerably  higher  than  what  is 
presently  available. 

In  any  event,  the  practicality  of  this  device  will  depend  upon  the 
reading  speed,  which  depends  on  the  exact  system  used  and  the  method  of 
reading. 

The  second  possible  use  of  electron  stimulated  dissociation  is  in 

the  growth  of  a material  on  a substrate,  e.g.,  silicon  on  silicon.  In 

this  case  a layer  of  silan^  (SiH^)  would  be  continually  deposited  onto  a 

silicon  surface  followed  by  local  irradiation.  Upon  irradiation  the  silane 

•/ 

dissociates  leaving  a silicon  outcropping.  By  continually  replacing  the 
silane  layer  a local  deposit  is  formed.  Using  very  small  diameter 
electron  beams,  deposits  of  small  dimensions  could  be  made.  A process 
such  as  this  could  be  used  to  fabricate  surface  gratings  and  integrated 
circuits  of  a size  smaller  than  that  which  is  currently  possible. 
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APPENDIX  B 
Conputer  Programs 

This  appendix  is  composed  of  four  computer  programs.  The  programs 
are  all  written  in  FORTPAN,  and  each  has  a brief  description  at  the 
start . 


I 
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LEED.F4 


Program  LEED.F4  is  used  to  calculate  the  interference 
function  for  the  data  listed  in  file  G00F1.DAT,  alon?  a 
characteristic  line.  That  line  is  determined  by  setting  AH 
(AH=(b/X  )cos  3)  equal  to  a fixed  value  and  varing  AH2 
(AH2=(a/\  )cos  a)  about  the  area  of  interest. 
DIMENSION  I(32,^2) ,F(32,32) 

CALL  IFILE(20,  'G00F1  ' ) 

90  READ(20,100)  I 

100  F0RMAT(32I) 

130  PI2=6 .2831753071^96 

CALL  0FILE(20,  'DEP1  ' ) 

N = 240 

WRITE  (20, 1 1 )iM 

11  FORMAT(I) 

DO  12  JJ  = 1 , M 
AJJ=JJ 

AH=AJJ*.01-1  . 15 
WRITE(20 , 13)AH 
13  FORMAT(F16.8) 

12  CONTINUE 

DO  270  JJ=1  ,N 
AJJ=JJ 

45  AH  = AJJ*.01-1  . 15 

46  AH2=0.0 
140  SUMIrO. 

145  SUM2=0. 

150  DO  240  M=1 ,32 
160  DO  230  N=1 ,32 

AG=PI2»( (N-16)*AH+(M-16)*AH2) 

F ( M , N ) = I ( M , M ) 

XX=F(M,N) 

190  A=(XX)»C0S(AG) 

200  E=XX»SIN(AG) 

210  SUVUSUMUA 
220  3Uy2=3U'^2  + 5 

230  CONTINUE 
240  CONTINUE 

245  3UM=3UM1 **2.+3UM2»»2. 

IF(AJJ.GE.2)  GO  TO  1233 
WRITE(20,  123)N 

123  FORNATCG) 

1233  WRITF(20, 124)SUM 

124  FORMAT(G) 

270  CONTINUE 

290  END 
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LDIF.F4 


This  DroiJra.T  calculates  the  variation  of  the  Aueer 
intensity  across  a line  of  width  2H  in  accordance  with 
Eo.  (51),  for  a fixed  value  of  Dt.  This  calculated 
intensity  is  then  comoared  to  the  exoeriTental  obtained 
values  of  '^(t)/Mo  and  Dt  is  changed  until  the  best 
least  souares  fit  is  obtained,  from  which  the  diffusion 
coefficient  D is  obtained.  The  incut  variables  are  the 
Auger  oeak-to-peak  height  PPO(IK)  and  the  position 
XX0(I!<)2.  The  output  (Dt)  is  AMIN. 

DIMENSION  A(25)  ,B(30)  ,EP(4)  ,C(4) 

DIMENSION  ST(1000),XX1(13),PP0(1?) 

CALL  IFILE(20,  'LC851  ’ ) 

NNr 1000 


22 

24 


23 

27 


N = 9 

DO  24  IK  = 1 ,N 
READ(20,22)XX1 (IX) 
F0RMAT(F14. 6 ) 
WRITE(5,22)XX1 (IK) 
CONTINUE 
DO  27  IK=1 ,N 
READ(20,23)PP0(IK) 
F0RMAT(F14.8) 
VRITE(5,23)PP0(IK) 
CONTINUE 


SSDT=2.43E-2 
v = 24 

A(1 )=16443152242714E-13 
A(2)=-49760497548E-13 
A(3)=  0643570883797E-13 
A(4)=  0196418177368E-13 
A(5)=-0001244215694E-13 
A(6)=-C009101941905E-13 
A(7)=-Q001796219835E-13 
A(8)=  0000139836786E-13 
A(9)=  164789417E-13 

A(10)=  390C9267E-13 

A(n)  = - 008931  45E-13 

A(12)=-  03747896E-13 

A(13)=-  01298818E-13 

A(14)=-00136773E-13 
A(15)=  00077107E-13 

A(16)=  00046810E-13 

A(17)=  00011844r-i3 

A(18)=-  0005E-13 

A(19)=-  1384E-13 

A(20)=-  0652E-13 

A(21)=  0145E-13 

A(22)=  0010E-13 

A(23)=  0024E-13 
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A(24)=  0011E-13 

A(25)=  0002E-13 

130  CO  4000  IDT=1 , 1000 

SQDT=SODT+( re-6 )*IDT 
20  H=3. *.05/3.75 

DO  1000  IK=1,N 

30  X1=(2.67/2. 1 )»(XX1(IK)+(15.E-3) ) 

60  X1=X1-H 

70  X2=X1+2.»H 

0(1 )=(H-X2)/(2. »SODT) 

C(2)=(H-X1 )/(2. *SQDT) 

C(3)=(H+X2)/(2. *S0DT) 

C(4)=(H+X1 )/(2. *S0DT) 

SOIPIz. 56418958 
E=SODT/(2. »H) 

IJ=0 

80  IJ=IJ+1 
90  W=C(IJ) 

X=A3S  (V) 

IF  (ABS  (X)-.01)  1,  2,  2 

1 XERR=2. 0/(3. 0*1 . 77245385 )*X«( 3. 0-X»*2) 

GO  TO  6 

2 Z=(X-1 .0)/(X+1 .0) 

DO  3 1=1,30 
E(I)=0. 

3 CONTINUE 
DO  4 1 = 1 ,M 
Ml =(M+1 )-I 

R(M1  )=2.0*Z»B(M1  + 1 )-B(mi+2)+A(MU1  ) 

4 CONTINUE 

F = -E(2)  + Z»P(1 ) + .5»A(1  ) 

XEPRrI .0-(  1 .0/1 .77245385  )»( EXP  (-(X*»2)  ) )»F 
IF  (ABS  (X)-.01)  6,  7,  7 

6 CERP=1 .0-XERR 
GO  TO  5 

7 CERP=(1.0/1.77245385)»(EX?  (-(X»»2)))»F 

5 ERR=XERR 

JI=2 

IF(W)200,210,210 
200  JI=1 

210  ER(IJ)=ERR»(-1 )»»JI 
IF(IJ-4)80, 150, 150 

150  AMTT  = -C(1  )»ER(1  )+C(2)*ER(2)-t-C(3)*ER(3)-C(4)» 

ER(4) 

TT=SQIPI«(EXP(-C(2)»C(2) )+EXP(-C(3)*C(3) )-E 
XP(-C(1 )*C(1 ) ) ) 

TTT=S0IPI*(-EXP(-C(4)«C(4) )) 
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AMT=E»(AyTT+TT+TTT) 

170  STA=( (AVT-PPOCIK) )**2. ) / ( ( -9 . 857E-2 ) » . 75 »PPO 
(IK)+9.957E-2)»*2. 

ST(IDT)  = STaDT)+STA 
CHI2=ST(IDT) 

1000  CONTINUE 

77  INK=IDT-1 

IFCIDT.EO. 1 ) GO  TO  4000 
AyiNL  = SCRT(ST(INf<:)  ) 

IF(ST(IDT).LT.ST(IN’-<))GO  TO  100 
GO  TO  3000 

100  AyiN2=30RT(ST(IDT) ) 

CO=ink 

4000  continue 

3000  DDD=SSDT+(1E-6)»(CQ+1 . ) 

WRITE(5,560)AyiN2,DDD 

560  FCRyAT(7H  AmIN2= , E 1 6 , 3 , 1 OH  SORT ( DT  ) = , E 1 6 . 8 ) 

5000  END 


1. 
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P00.F4 


This  oroerar"  is  used  to  calculate  the  theoretical  Aueer 
intensity,  in  accordance  with  Ec.  (53)  when  rraterial 
is  lost  by  both  diffusion  and  desorption  f rorri  a line 
source.  The  values  of  the  activation  energy  of 
desorption,  QDES , diffusion  ,ODIF  and  temoerature  (T5) 
are  varied  until  the  best  fit  is  obtained. 

DIMENSION  C(0/50Q0)  ,0(0/5000) 

PEAL  LI ,K,M2 
Z=2»*12 
T = 0 
T1  =0 

0DIF  = 215iJ7. 

QDES=80000. 

T5=1000. 

D1  = (ii.66E-2)»EXP(-ODIF/(  1 .9872«T5)) 

X2= 1E-3«/D1 

:<  = (il.6/39.6)»EXP(-ODE.S/(  1 .9872*T5)  ) 

XI =SQRT( 1E-34/K) 

C1  = 1 . 188E1U 
C0=2.772E14 
C2=C1+C0 
S = 5 

LI =.05625 
T0=1 
H1=6 
H2  = 3 
H = 10 

T9=1000. 

CALL  0FILE(20 , 'P0U8  ' ) 

K1=2*»H1 

K5=K1+1 

K2=2«»S 

P=K2»3. 1^15926 
H3=(L1/X1  )**2 
H4=T0/2*»H2 
H5=1/H 

WRITE(20, 10) 

WRITE(20, 15)  ,D1  ,t(,S,H1 
WRITE(20,20) 

t^RITE(20,25)  ,C0,C1  ,H2,L1 
VRITE(20,26) 

WRITE(20,27) ,0DIF,0DES ,T5 


10 

F0RMAT( ' D K S 

HI  ' ) 

15 

F0RMAT(4G ) 

20 

FORMAT ('  CO  Cl 

H2  LI ' ) 
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25 

26 
27 


100 


120 


1U0 

30 

C(3L/2)  ' ) 
500 


532 

160 

5^1 

180 


6 1 6 
612 

619 

622 

200 

631 

220 


FORMAKiiO ) 

formak'  odif  odes  temp,') 
F0RMAT(3G  ) 

DO  100  M = 0,:<1 

C(N)=C2 

CONTINUE 

C(N)=C1+CO«(1+COS((P«(N-K1  ) )/FC1  ) )/2 

CONTINUE 

DO  140  fJ  = K9,Z 

C(N)=C1 

D(N)=0 

CONTINUE 

WRITE(20,30) 

FORMATC'  TIME  C(0)  C(L/4)  C(L/2) 

D(0)  = 2»(C( 1 )-C(0)  )/H3 
DO  160  N=1 ,Z 

DCN)  = (C(N+1 )-2»CCN)+C(N-1  ) )/H3 
IF  ((C(M-8)-C(N)).GE. IE-6)  GO  TO  160 
IF(N.GE.K5)  GO  TO  532 
GO  TO  160 
J = N 

GO  TO  541 
CONTINUE 
DO  180  N=J,Z 
C(N)=C(N-8) 

CONTINUE 
DO  200  N=0,Z 
A5=C(N)»K 

IF(ABS(C(N) ) .LT.X1  ) GO  TO  612 
IF(ABS(D(N)) .LT.X2)  GO  TO  616 
C(N)=C(N)+H4*(D1 »D(N)-A5»C(N)  ) 

GO  TO  619 

IF(ABS(D(N)) .LT.X2)  GO  TO  619 
C(M)=C(N)+H4*D1 *D(N) 

GO  TO  619 

C(N)=C(N)-H4*a5*C(N) 

IF(C(N-8)-C(N)  ,GE. IE-6)  GO  TO  200 
IF(M.GE.K5)  GO  TO  622 
GO  TO  200 
J = N 

GO  TO  631 
CONTINUE 
DO  220  N=J,Z 
C(N)=C(N-8) 

CONTINUE 

T=T+H4 


C(L) 
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T1=T1+H4 

IF(T1 .LT.H5»T9)  GO  TO  500 

WRITE(2n,50)  ,T,C(0)  ,C(  .5*Kl/2.  ) ,C(K1  /2.  ) ,Cf:< 
1 ) ,0(3.  *!C1/2.  ) 

50  FORMAT(6E12.6) 

T1=0 

IF(T.LT.T9)  GO  TO  500 
WRITE(20,60) 

60  FORMATC ' X AREA  ’ ) 

DO  260  1=1,30 
M1=(I-1 )»K1/8. 

E = 0. 


I — w , 

A2  = 0. 

EU  = 0 . 

DO  240  M = N1  + 1 ,M1+f(1 
A = C(N)/(K1*2*(C0+C1  )) 

E = E+A 

240  CONTINUE 

IF(N1 .GT.K1-1  ) GO  TO  870 
GO  TO  920 

870  DO  280  N=N1-K1 ,N1 

A = C(N)/(2*K1 »(C0  + C1  ) ) 
FPrFP^-A 

280  CONTINUE 

GO  TO  1000 

920  DO  300  N=0,N1 

A = C(N)/(2*fC1  »(C0  + C1  ) ) 
E3=E3+A 

300  CONTINUE 

DO  320  M=1 ,K1-N1 

A = C(M)/(2*K1 »(C0  + C1  ) ) 

E4=E4+A 

320  CONTINUE 

1000  A1=E+E2+E3+E4 

N2  = N1*L1 /(K1  ) 
VRITE(20,55),N2,A1 
55  FORMAT(2G) 

260  CONTINUE 

END 


150 


FLADES.F'^ 


THIS  PROGPAw  SEPERATES  THE  FLASH  DESORPTION  SPECTRUM 
INTO  THE  DESORPTION  OF  TWO  FIRST  ORDER  DESORPTIONS . THE  BEST 
ESTIMATE  OF  THE  FIRST  Tp  IS  PUT  IN,  THE  OUTPUT 
IS  El  OF  THAT  PEAK. THE  INPUT  DATA  IS  N/Mn  AND  TEMP 
AT  THE  VALUES  OF  N/Np.THE  OUTPUT  IS  THE  NORMALIZATION 
OF  THE  SECOND  CURVE, ITS  TpCGIVEN  ASTP2)  AND  E2. 

DIMENSION  T( 12)  ,Y( 12) ,YTHE ( 1 2 ) , Y 1 THE ( 1 2 ) ,Y2THE ( 12) ,TE 

MP1 (12) 

DIMENSION  TFMP2(12) , CHISO ( 200C0 ) , Y 1( 1 2 ) , Y2 ( 1 2 ) ,CHISOO 

(20000) 

10  CALL  IFILE(20,  ’DES2  ' ) 

15  DO  27  1=1,9 
20  READ(20,25)  Y(I) 

25  FORMAT(G) 

27  CONTINUE 

28  DO  37  1=1,9 

30  PEAD(20,35)  T(I) 

35  FORMAT (G) 

T(I )=T(I )+20. 

37  CONTINUE 
TP=172. 

DO  Hi!  1 = 1,5 

WRITE(5,^3)  Y(I) ,T(I) ,TP 
«3  FORMATdH  ’N/NP=  ' ,G  , 'TEMP=  ' ,G  , 'TP=  ' ,G) 

Uii  CONTINUE 
50  R=1.9g72 

80  E0=8^75. 

90  E00=7800. 

1 10  DO  250  11=1,20 
125  CHSO=0.0 
120  DO  210  1=1,5 
130  E1=E0+50.»II 

133  X=E1/(TP*R) 

135  XX=E1/(T(I)»R) 

1U0  TEMPI (I )=XX-X 

150  YT=TEMP1 (I)-1 . + (EXP (-TEMPI ( I ) ) ) * ( (T ( I )/TP ) »»2 . ) 

160  YTHE(I)=1 ./(EXP(YT) ) 

170  CHI=(Y(I)-YTHE(I))*«2. 

180  CHSO=CHI+CHSO 
190  IF(I.LE.a)  GO  TO  210 
200  CHISOOdI  )=CHSQ 
210  CONTINUE 
220  IFdI.EQ.  1 ) GO  TO  250 
IA=II-1 
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2^0  IF(CHISOC(IA).GT.CHISCO(II))  GO  TO  250 
2U0  En=F0+50.  *(IA) 

ABC  = CHISCQ (lA  ) 

GO  TO  260 
250  CO^'TI^fUE 
260  WRITE(5 ,270)  Ell ,APC 

270  FORMATdH  ’E1  = ',G,'  CHISO=’,G) 

IC0UMT=0 
290  T0=185. 

IJKArO 
CHSO=0.0 
ROO  IJK=0 
310  DO  570  JJ=1 ,U0 
320 • DO  560  JJK  = 1 , 10 
330  DO  550  JJM=1,20 
CHSQ=0.0 

340  DO  490  J=3,9 
350  TP2=T0+JJK 
360  .E2  = E00+50.  ♦JJ 

370  X = E2/(R»T(J)  ) 

380  XX=E2/(R»TP2) 

390  TEMP2(J)=X-XX 

XY  = E1 1/(R*TCJ)  ) 

XXY=E1 1/(R»TP) 

TEMPI (J)=XY-XXY 

400  Y1TH  = TEMP1  (J)-1  . ( EX  P ( -TEW  P 1 (J)))*((T(J)/TP)««2.) 

Y1THE(J)=1 ./(EXP(YITH)) 

410  Y2CJ)  = (Y(J)-(Y1THE(J)  ) ) 

420  Y2TH=TEMP2(J)-1  . + CEX?  ( -TE‘'P2  ( J ) ) ) » ( CT  ( J ) /TP2  ) « *2  . ) 
430  Y2THE(J)  = 1 ./(EXP(Y2TH)  ) 

4 40  CHI  = (Y2(J)-(Y2THECJ) )»( .R8f JJM».01  ) )**2. 

450  CHSQ=CKI+CHSC 

460  IFCJ.LE.B)  GO  TO  49c 

470  IJK=IJF+1 

480  CHISQ(IJK)=CHSO 

490  CONTINUE 

500  IFCJJN.EO.  1 ) GO  TO  550 
IIN=IJK-1 

510  IF(CHISO(IIM)  .GT.CHISO(IJK)  ) GO  TO  550 
IFCIJKA.LT. 1 ) GO  TO  520 

515  IFCCHISQdJKA)  .GT.CHISOCIIN)  ) GO  TO  520 
GO  TO  550 

520  AJJN  = ( (JJN-1  )».0U.38  ) 

530  AJJK=JJK+TO 
540  AJJ=JJ»50.+E00 
IJXArllN 
GO  TO  560 
550  CONTINUE 
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560  COMTIMUE 
570  CONTINUE 

580  WPITE(5,590)AJJN ,AJJK  ,AJJ 

590  FORNATdH  ' NORN  . CONST  . = 0 TP2  = ’,  G , ' EP2=',G) 

670  IJKArO 
680  IJKrO 
INKrO 
CHS0  = 0 .0 

690  DO  1070  JJ  = 1 ,iJO 
700  DO  1060  JJK=  1,10 
710  DO  1050  JJN=1 ,20 
600  IJ=3+INK 
IF(IJ.GT.9)  GO  TO  780 
610  X=(AJJ/R)*(1 ./T(IJ)) 

XX  = (AJJ/(R»AJJK)  ) 

TEMP2(IJ )=X-XX 

620  Y2TH=TEMP2(IJ)-1 . + (EXP ( -TEMP2 ( TJ )))«(( T ( IJ  ) /A JJK  )« »2 . ) 
630  Y2THE(IJ)=(1 ./(EXP(Y2TH)))»AJJN 
6U0  WPrTE(5,650)  Y2THE ( IJ  ) , T ( IJ  ) 

650  FOPMATdH  'N/NP=',G,’  TEMP=',G) 

INK=INK+1 

IF(IJ.E0.9)  GO  TO  780 

GO  TO  600 

780  CHSO=0.0 

790  DO  930  J=1 ,6+ICOUNT 

800  TP1 =TP+JJK-P . 

810  E12=E1  1-1006.+25dJJ 
820  X = E12/(R»T( J)  ) 

830  XX=(E12/(R*TP1  ) ) 

8^0  TEMPI (J)=X-XX 

850  Y1TH=TEMP1( J)-1 . + (EXP ( -TEMPI ( J ) ) ) * ( ( T ( J ) /TP  ) * *2  . ) 

860  Y1THE(J)=1./EXP(Y1TH) 

870  Y1 (J)  = (Y(J)-Y2THE(J) ) 

880  CHI=(Y1(J)-(Y1THE(J))»(.82+.01»JJN))»*2. 

890  CHSOrCHI-^CHSQ 

900  IF(J.LE.5+IC0UNT)  GO  TO  930 

910  IJK=rJK+1 

920  CHISO(IJK)=CHSO 

930  CONTINUE 

IFdJK.EQd  ) GO  TO  1050 
950  IINrlJK-l 

960  IF(CHISO(IIN)  .GT.CHISO(IJK) ) GO  TO  1050 
970  IFCIJKA.LT. IIGO  TO  1000 

980  IFCCHISOdJKA)  .GT.CHISO(IIN)  ) GO  TO  1000 

990  GO  TO  1050 

1000  A JJN=(JJN-1)*.01+.82 

1010  AJJK=JJK+TP-3. 

1020  AJJ  = JJ»25. ■►El  1-1000. 

1030  IJKArllN 
lOiiO  GO  TO  1060 
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1050  CONTINUE 
1060  CONTINUE 
1070  CONTINUE 

1080  WRITE(5, 1090)AJJN,AJJK,AJJ 

1090  FOPVATdH  'NORM. CONST  1 = ’ , G , ' TP  = ' , G , ' EP1  = ',G) 

1100  DO  1180  J=1 ,6+ICOUMT 
1110  X=(AJJ/R)»(1 ./T(J)) 

1120  XXr  (AJJ/(R*AJJf()  ) 

1 130  TEMPI (J)=X-XX 

11«0  Y1TH  = TEMP1(J)-1 . + (EXP(-TEMP1 (J)) )»((T(J)/AJJK)»*2.  ) 
1150  Y1THE(J)=(1 ./(EXP(Y1TH)))*AJJN 

1160  writE(5, 1 170)Y1THE(J)  ,T^J) 

1170  FORMATCIH  'N/NP1=',G,'  TEMP=’,G) 

1180  CONTINUE 
ICUUNT=IC0UNT+1 
IF(IC0UNT.E0.3)  GO  TO  1190 
El 1=AJJ 
GO  TO  290 
1190  END 
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APPEND  LX  C 
Error  Analysis 


The  error  analysis  of  the  various  experimental  parameters  is  straight- 
forward. In  our  experiments  we  generally  want  to  determine  the  quantity 
X = f (u,v,  ...)  by  measuring  u,  v,  ...  . Although  it  may  not  be  exact, 
we  assume  that  the  most  probable  value  of  x is  given  by  [45], 


X = f (u,v,  . . .)  , (C.l) 

The  uncertainty  in  the  resultant  value  of  x can  be  found  by  considering 

the  spread  in  x resulting  from  combining  the  individual  measurements 

u.,  v.,  ...  into  individual  results  y . , 
i’  1 i 

Xf  = f (u^,v^,  . . .)  . (C.2) 


The  variance  of  x,  ^ , 

X 


is  given  by. 


X 


N 

lim  i 2 (x  -x)^ 
N-»  ^ i=l  ^ 


(C.3) 


The  deviation  x.-x  can  now  be  expressed  in  terms  of  the  deviations 
1 

- u,  v^  - V,  or, 


X^  - X 


/ dx  , . dx  . 

('^i  ■ d[r  + d7 


(C.4) 


Combining  Eqs . (C.3)  and  (C.4)  yields 


C ^ lira  - r ; (u.  -u)  3—  / + (v.  -v)  ‘ — I + . . . ! 

X „ N . - 1 1 Vdv'  j 

N-  * 1 


lira  ^ eI(u.-u)V^)  + (v.-v)^(^) 

„ ^ N . _ 1 ^du/  1 Vdv/ 

N -*  ® 1 


(C.5) 


2 — 2 2 — 2 
Notice  that  c =lim-Z(u.  -u)  and  7 =lim~E(v.  -v)  ; similarly 

“ s-=  “ 1 " ''  s-=.  1 ‘ 

2 1 =*  - - 

= lira  — E [ (u^-u) (v^-v)] . Equation  (C.5)  can  now  be  written  as 
N “ i 

^ (7^)  +cr^  (~)  + 2C^  (j^)  + ...  . (C.6) 

X u Vdu/  V '-dv/  uv  vdu/  vdv' 

If  the  fluctuations  of  u,  v,  ...  are  uncorrelated,  Eq . (C.6)  reduces  to 


2 _ 2 /'dxV  ^ 2 /dx)^  ^ 


(C.7) 


or  written  more  generally, 


2/‘'^jN 


(C.8) 


In  order  to  be  able  to  use  Eq . (C.8)  to  estimate  the  standard  deviations 

of  the  diffusion  and  desorption  parameters,  we  need  to  know  the  variance 

2 

of  each  measurement  cr. . When  using  Auger  spectroscopy  to  study  the 

2 

diffusion  and  desorption  of  nitrogen,  is  the  variance  of  the  nitrogen 
Auger  peak. 
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The  value  of  is  determined  by  first  depositing  a predetermined 

amount  of  electron  dissociated  nitrogen  on  the  surface.  A number  of 

Auger  scans  are  then  taken  in  the  region  of  nitrogen  K L,  , L„  . 

transition.  From  the  variation  in  the  peak  to  peak  height  of  the  nitrogen 

2 

Auger  transition,  a value  of  is  obtained.  We  found  that  can  be 

described  by  o.  = -7.39  X lO'^  PPO(IK)  = 9.957  X lO"^  where  PPO(IK)  is 

the  peak  to  peak  height  of  the  normalized  nitrogen  Auger  transition. 

As  an  example,  consider  the  estimation  of  the  standard  deviation  of 

the  diffusion  coefficient  D of  the  electron  dissociated  nitrogen  after 

heating  for  3500  seconds  at  800°K.  The  computer  program  LDIF.F4  is  run 

with  the  appropriate  values  of  XXl  (position)  and  PPO  (.Auger  peak  height). 

•8  2 

A value  of  D = 7.55  X 10  cm'^/sec  is  obtained.  One  of  the  values  of  PPO 
is  now  changed  by  a small  amount  (~  17o)  and  another  value  of  D is  calculated. 
This  procedure  is  repeated  for  all  the  data  points.  We  now  have  values 

dDi 

of  , and  Eq.  (C.8)  becomes 

d(PPO)^ 

2 r-  2 do 

a = I o.  , 

° i ^ d(PPO)^ 


from  which  the  standard  deviation  of  D can  be  obtained. 


Once  the  values  of  are  obtained  for  various  temperatures  they 


are  plotted  as  ln(D)  vs.  1/T  using  the  program  LEASTS. BAS.  This  program 
computes  the  best  least  squares  fit  to  the  data  with  a weighting  function 
w^  given  by  [46], 


‘JU. 

1 


(C.9) 
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I 
I 

An  identical  procedure  is  followed  for  the  calculation  of  the  errors 

I 

i I associated  with  the  desorption  parameters. 
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